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Gottesman, Kitaev and Preskill have formulated a way of encoding a qubit into an oscillator such 
that the qubit is protected against small shifts (translations) in phase space. The idea underlying 
this encoding is that error processes of low rate can be expanded into small shift errors. The qubit 
space is defined as an eigenspace of two mutually commuting displacement operators S p and S q which 
act as large shifts/translations in phase space. We propose and analyze the approximate creation 
of these qubit states by coupling the oscillator to a sequence of ancilla qubits. This preparation of 
the states uses the idea of phase estimation where the phase of the displacement operator, say S p , 
is approximately determined. We consider several possible forms of phase estimation. We analyze 
the performance of repeated and adapative phase estimation as the simplest and experimentally 
most viable schemes given a realistic upper-limit on the number of photons in the oscillator. We 
propose a detailed physical implementation of this protocol using the dispersive coupling between a 
transmon ancilla qubit and a cavity mode in circuit-QED. We provide an estimate that in a current 
experimental set-up one can prepare a good code state from a squeezed vacuum state using 8 rounds 
of adapative phase estimation, lasting in total about 4/r sec., with 94% (heralded) chance of success. 


I. INTRODUCTION 

It is to be expected that quantum error correction 
will be essential in the implementation of reliable, 
large-scale, quantum computation. Efforts are underway 
for superconducting transmon qubits jlj coupled to 
microwave-resonators (circuit-QED) to realize compo¬ 
nents of a surface code architecture in which robust 
logical qubits are comprised of many 0(10 — 100) 
elementary transmon qubits and resonators @, II- As 
the qubit and resonator overhead of such encoding is 
rather daunting and inefficient, one may ask whether 
there exist more efficient ways of using high-coherence 
qubits and long life-time microwave resonators (see 
some numbers in Table IU) to encode robust logical 
qubits. In this paper we will explore the proposal by 
Gottesman, Kitaev and Preskill [J] to encode a qubit 
into an oscillator and take a first stab at analyzing 
how a qubit can be encoded and preserved in a single 
mode of a microwave cavity. This scheme may not 
only be interesting as an alternative route towards 
scalable quantum computation, but also as a way of 
generating a highly non-classical state in a (microwave) 
cavity which one can actively protect against photon loss. 

The GKP (Gottesman-Kitaev-Preskill) code is of in¬ 
terest as it offers the possibility to use the cavity mode 
to store a single qubit, while a transmon qubit coupled to 
this cavity mode can be used to manipulate the state of 
the cavity. The goal is to prepare and preserve a qubit in 
a cavity mode such that it has a coherence time much 
longer than that of the cavity or the transmon qubit 
itself. The GKP scheme is attractive in that gates on 
this qubit are relatively straightforward as Clifford gates 
can be realized using linear optical componenents, see Q 
and Appendix [A] In addition, this qubit-into-oscillator 
scheme can be concatenated with a surface code or an 
other stabilizer code for scalable protection (see e.g. @ 


for a surface code with oscillators). 

However, a realistic scheme for preparing the encoded 
states and performing quantum error correction has been 
missing so far (in Section II Bl we list some of the previ¬ 
ous work) and this is the focus of the present paper. In 
the remainder of this section we review and discuss the 
code states and their highly-nonclassical properties. We 
argue why, at a heuristic level, using phase estimation 
can provide a means for approximately generating such 
quantum states. In Section E3 we review and introduce 
some tools and formalism which are useful in assessing 
the quality of a preparation protocol, using a model of 
displacement or shift errors. In Section fl Dl we physically 
motivate the shift error model by discussing how phys¬ 
ical errors can be expanded into such shifts. The shift 
error model plays a similar role as a Pauli error model 
for qubits. This section partially reviews some of the in¬ 
sights of 3; it adds to these by showing in Section IIP II 
how an expansion into shift errors can lead to a reduction 
in error rate by several orders of magnitude. Section HT1 is 
devoted to exploring and understanding how to prepare 
an approximate code state using the protocol of phase 
estimation, in this case determining the eigenvalue of a 
unitary displacement operator. There are many ways to 
do phase estimation leading to different approximations 
for the code states and using different numbers of pho¬ 
tons. These phase estimation protocols all consist of a 
coupling a sequence of single qubit ancillas to an oscil¬ 
lator mode sequentially and measuring the ancilla qubit. 
We argue why we focus on two simplest protocols: one 
is repeated phase estimation without feedback (Section 
III Bl) and a second one is a simple form of adaptive phase 
estimation with feedback (Section III Cl) . For those two 
schemes we explicitly show in Section III El how well with 
a small number of rounds/ancilla qubits an approximate 
code state can be prepared. We also consider the average 
number of photons used and the variance thereof, as these 
numbers are relevant in any experimental realization of 
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the protocol. 

In Section [TTT] we discuss a realization of the code state 
preparation protocol in an experimental circuit-QED set¬ 
up where the oscillator corresponds to a cavity mode. We 
review the circuit-QED transmon-cavity coupling and the 
physical strength of various parameters as reported in re¬ 
cent experiments in circuit-QED. We discuss several as¬ 
pects of an implementation such as the use of a tunable 
coupling and the implementation of a (transmon) qubit 
controlled-displacement gate by means of a controlled- 
rotation or a direct pulse-driven controlled-displacement 
(details in Appendix D). In a final Section fill Al we dis¬ 
cuss at a qualitative level the dominant sources of inac¬ 
curacy in implementing the protocol; we especially note 
that the presence of cavity nonlinearities may require a 
deeper analysis. We end the paper with a discussion and 
outlook for future work. 

In Appendix A and B we collect some background tech¬ 
nical results which are relevant for implementing gates on 
the code states and performing quantum error correction. 
Appendix C gives background on the choice of feedback 
phases in the phase estimation protocol in Section III Cl 


A. GKP Code States 

In order to introduce the code states, we assume a har¬ 
monic oscillator with annihilation operator (creation op¬ 
erator) a (a') with which one can define (dimensionless) 
quadrature operators p = ^(aJ — a) (momentum) and 

q = -^(a + a^) (position) obeying the canonical com¬ 
mutation relations [q,p\ = %■ The ideal GKP code is 
simply a prescription of a two-dimensional subspace in 
the infinite-dimensional oscillator space. This subspace, 
called the code space, is defined as the +1 eigenspace 
of two commuting displacement operators S p = e ~ l2 ^ p 
and S q = e l2 '^ q (One can verify the commutation of 
these operators by using the identity e A a B = e B e A e^ A,B ' i 
for A , B linear combinations of p and q). The space fixed 
by S p and S q is two-dimensional as there are additional 
operators which we can identify as X and Z which com¬ 
mute with both S p and S q , but which mutually anti- 
commute. 

States in this two-dimensional subspace will be called 
code states as they encode a qubit. Of course, there 
are many possible ways of choosing a two-dimensional 
subspace of a harmonic oscillator to define a qubit, e.g. 
choose two energy eigenlevels or pick two orthogonal cat 
states to define |0) and |1). The subspace of the GKP 
code is one in which |0) and |1) are highly non-classical 
states and is chosen such that small phase space displace¬ 
ments can be undone, by quantum error correction. In 
this sense the code offers the possibility to realize a long- 
lived, well-protected qubit. 

What is immediately interesting about a state in this 
subspace is that the condition S p = 1 and S q = 1 fixes 
both momentum p = 0 mod y/n and position q = 0 


mod y/n to be sharply determined. A code state thus 
escapes Heisenberg’s uncertainty principle by being not 
localized at a single p and q but being a superposition 
of many equally spaced sharp values for p and q. Due 
to having low variance in both quadratures, we can ex¬ 
pect that a code state is useful in metrology for detecting 
small displacements which shift p or q by less than y/n. 

If we write the operators S p and S q as displacements, 
one has S p = D(y/2n) and S q = D(iy/2n) (using D(a) = 
exp(aa^ — a*a) with D{a) |vac) = |a) and |a) a coherent 
state). In this code space one can define the qubit states 
|0) and |1) as 

oo oo 

|°) oc Y D (tVEr) \q = o) = Y s p \ q = °) ’ 

t=— oo t=— oo 

oo oo 

|i)oc Y D{tVEH:)\q = y/n) = Y St p\q = V^)- 

t=— oo t=— oo 

(1) 

Thus, |0) is a uniform superposition of eigenstates of 
position q at even integer multiples of y/n, while |1) 
is a uniform superposition of eigenstates of position q 
at odd integer multiples of y/n. One can also consider 
these states as superpositions of p-eigenstates, i.e. |0) oc 
I d P (EZ-oo S t ) \p) oc EZ-ao I P = S V^) while |1) oc 

I dp e’^Et.oo 3 p)Ip> « E^L-ooM)" \p= S ZZ 

(using that EZ-oo e~ l2 ^ pt oc E7=-oo S (p = 

Thus, both |0) and |1) have nonzero amplitude at integer 
multiples of y/n inp-space, but due to alternating phases, 
these amplitudes destructively interfere at odd multiples 
of y/n in the state |+) = ^(|0) + |1)) [f|. The opera¬ 
tor X: |0) |1) is given by X = e~ l ^ p , as it shifts 

each eigenstate |<?) by y/n. Note that X 2 = S p which 
equals I only on the code space (+1 eigenspace of S p ) 
(see discussion in Appendix 0. One can write 

oo 

|+} = -=(|0) + |l»cx Y S \\P= °) 

^ t= — OO 

OO 

|-) = - 7 =(|0)-|l))oc Y S l\ P=V^>, (2) 

t=— OO 

with Z: |+) o |-) given by Z = e^ q and Z 2 = S q . 

A simple way to understand the ideal preparation of 
a code state |0) (or |1}) is through Eq. (JT]) : one starts 
with a +1 eigenstate of Z, namely \q = 0) to which one 
applies n Sp=1 = Efe-ocL^p which is the projector onto 
the +1 eigenspace of S p p}. 

1. Approximate GKP Code States 

Naturally, Ref. Q realized that the perfect code states 
in Eq. 0 are unphysical as their preparation would re¬ 
quire infinite squeezing and an infinitely sharp projection 
onto the +1 eigenspace of S p . Ref. Q suggested using 
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approximate code states: we can understand how this 
approximation comes about in, say, the definition of |0), 
as follows. 

As it would take an infinite amount of squeezing to pre¬ 
pare \q = 0) we replace this state by a finitely squeezed 
state (in q ) centered around q = 0, that is, a squeezed 
vacuum state, Eq. ©■ Then we need to implement an 
approximate projection onto the S p = 1 eigenspace. Such 
projection could come about by approximately estimat¬ 
ing the eigenvalue of the operator S p and then post- 
selecting the measurement outcome on this eigenvalue 
being +1 (while discarding other results). Such approxi¬ 
mate projection by post-selection, using an ancilla qubit, 
has been considered in Q, without making reference to 
any particular technology. However, any post-selection 
scheme will have a low probability of success and is in 
fact unnecessary. If one can implement a highly accurate 
phase estimation measurement of the unitary displace¬ 
ment operator S p , estimating its eigenvalue as some e l6 , 
then one can also correct for such an eigenvalue shift and 
shift, or displace back to the +1 eigenvalue code space. 

The formal definition of the approximate states is as 
follows. One replaces each delta-function in position q by 
a squeezed Gaussian state while the uniform superposi¬ 
tion over these localized states is replaced by a Gaussian 
envelope. For example, for the approximate |0) state, 
one starts with the squeezed vacuum state |sq.vac) with 
squeezing parameter A 2 = e~ 2r 

Isq.vac) = J ^-A_ e -?7(2A 2 ) | g) s ( 3 ) 

to which one applies a sum of displacements, as in Eq. CD. 
with a Gaussian filter, i.e. 

OO 

1°)approx E e~ 27rA2f2 D(tV 2^) |sq.vac) 

t=— OO 


A are identical, nor that the filter is Gaussian: in 0 a 
very general form of approximate code states depending 
on two different filter functions was formulated. 

One can obtain the form of the code states in p-space 
by taking a Fourier transform. For example, the state 

1+) = 72 (l°)approx + l^approx) haS the following form 

when A /y/n, -C 1, see Q: 

OO r. 

Ifo)approx rS E /e- AV/2 e-^- 2t ^) 2 /( 2A2 )|p)dp 

t=- OO ^ 

OO r. 

« E /e- 2 - A2t2 e- (p - 2 ^ )2 / (2A2) b)dp, 

t=- OO ^ 

(5) 

where one observes that the roles of A and A are inter¬ 
changed. If one encodes in a bosonic mode, it is natural 
to choose similar approximations in p and q as free evo¬ 
lution of the state evolves these quadratures into each 
other. This means that the choice A = A is natural. 
The filter is Gaussian, so that we can use squeezed states 
which are Gaussian in their quadrature spread. It can be 
seen in Eqs. © and © that the wave functions of ap¬ 
proximate code states are even in q and p , respectively. 
As the photon parity operator P = e l7rata transforms 
q —> —q, p —> —p, these code states will always have an 
even number of photons (this is not true for the approx¬ 
imate shifted code states in Eq. ©)• 

For an approximate code state with parameter A one 
can explicitly calculate the mean number of photons 
n = (ala) and its variance er 2 (n) = {(n — a^a) 2 ). Us¬ 
ing that n = {\{p 2 + q 2 ) + 1) = {q 2 ) + 1, A « 1/A and 
approximating the infinite sum over t by an integral, one 
obtains n rts —p- • A similar calculation yields 

v( n ) « ^ ~ fo ( 6 ) 


OO 

= Ej 

t=- OO J 

f e- 2 ^ 2 t 2 e-^- 2 t ^ 2 /^\q)dq. 

OO 

|1) approx 0(1 y! 1 

t —— OO J 

1 " e -nA 2 ( 2 t+l) 2 / 2 e -(q-( 2 t+l)^) 2 /( 2 A 2 ) 


(4) 


to leading order in 1/A. This shows that the approxi¬ 
mate code states are highly non-classical states with large 
fluctuations in photon number, i.e. scaling with the total 
| q) dqnumber of photons while the approximation parameter 
A controls the number of photons in the state. We will 
return to the number of photons in the approximately- 
prepared code states in Section Hi E II 


The code state is thus a Gaussian-weighted sum of dis¬ 
placed squeezed vacuum states. This code state and its 
Wigner function W(p,q) = 1 /7r e 2lpx d'*(q + x)d'(q — 
x)dx, where 'F(g) is the wave function of the state 
|0)approx * n th e position basis, are depicted in Figs. |T| 
and [2] for A = A = 0.2 which corresponds to 8.3 dB. 
Here and elsewhere in the paper we calculate dB as 
101og 10 (G) = 101og 10 (cosh 2 (r))) where G is the amount 
of gain through the amplifier and r = log(l/A) [g]. 

In the limit of A —> 0 (infinite squeezing) where 
the width of the Gaussian envelope becomes arbitrarily 
broad, one obtains the perfect code states. It is worth 
noting that it is not essential that the parameters A and 
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FIG. 1. Absolute value of the wave function of the approxi¬ 
mate code state |0) approx (a) and |l) approx (b) for A = 0.2. 
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FIG. 2. (Color online) Surface plot, contour plot and inte¬ 
gral over p , q of the Wigner function W ( p , q) for code states 
l°)approx with A = 0.5 ( a ) and A = °- 2 (b)- For q being 
an odd multiple of ypii, the Wigner function rapidly oscil¬ 
lates between positive and negative values when p is an even 
resp. odd multiple of y/jr. Integrating over the p-variable thus 
leaves, due to destructive interference, little amplitude for q 
at odd multiples of y/n, as can be seen in Fig. [T] 


B. Previous Work 

The preparation of the code states has been discussed 
in several papers see e.g. 0000. in 0 it was noted 
that an interaction H lnt = qb^b where q is the position 
operator of the oscillator (e.g. a cooled microinechanical 
system) in which the code state is to be encoded and b 
is the annihilation operator of a different bosonic mode, 
would be useful in preparing the code states. The pro¬ 
tocols discussed in |K. 12| use post-selection in preparing 
the code states making the success probability extremely 
small, while we find that (almost no) post-selection is re¬ 
quired by understanding the preparation step as phase 
estimation (see details in Section HlEl) . 

In 0 the use of cat states |0) oc |a) + |— a) and 
|1) oc | ia) + | —ia) was explored for encoding a qubit 
in a cavity mode and dissipative dynamics was consid¬ 
ered to preserve the even photon number code space. 
Even though dissipative dynamics can drive a state to 
this code space, it does not imply that such dissipative 
dynamics realizes quantum error correction as it has to 
approximately drive an encoded state in this subspace to 
itself. Unlike the approximate GKP code, the cat state 
code by itself (implemented in circuit-QED in 00) is 
not fully protected against loss of photons. The detection 
of a photon number change by 1 constitutes a reversible 
change (which can be undone), but when one does not 
detect this change, the cat states decay irreversibly, nor 
is the code protected against the simultaneous loss of two 
photons. 

In 0 multi-component cat states were dynamically 
generated in a cavity-mode using the Kerr nonlinearity 
of the cavity (see e.g. the description of such a protocol 
in 0). A multi-component cat state is a superposition 
of coherent states spaced at equal distances on a circle in 
phase space while the approximate code states are cre¬ 


ated by translations along a line. It seems hard to define 
two approximately orthogonal qubit states using these 
multi-component cat states such that (1) the phase space 
amplitude peaks for |0), |1) are sufficiently separated (so 
as to protect against small shift errors) and (2) |+) and 
|—) are similarly superpositions of sufficiently separated 
peaks. 


C. Quantum Error Correction and Protection 

The primary reason for defining the two-dimensional 
code space as the +1 eigenspace of S p and S q is that 
the code states are protected against small shift errors in 
phase space, of the form e _mp and e~ lvq with |w|,|u| < 
y/ir/2. Such shifts move the peaks of the code state in 
Fig. [Has e~ lup |g) = \q + u) and e~ lvq \p) = \p — v). Thus 
when the shift error is less than half the shift represented 
by the operators X or Z in magnitude, one can undo the 
error by shifting back the state by the minimal amount. 

Imagine that a shift error e~ lvq has occurred and the 
eigenvalue of the check operator S p is perfectly measured 
by some means. One will estimate the eigenvalue of S p 
as e l2 ^ v (as S p e~ mq |-0) = e l2 V* v e -™q |^,) for a state 
il> in the code subspace, as follows from the commuta¬ 
tion between e~ lvq and S p ). For |u| < y/n/Z, the phase 
8 = 2 v / 7ti> £ (—7r,7r) is uniquely given and thus one can 
correct the shift error by learning the phase 8 (and shift¬ 
ing back). For larger |u|, 8 is consistent with two values 
of v and choosing the wrong one leads to a logical qubit 
error. Error correction for shift errors e _mp works simi¬ 
larly by measuring the eigenvalue of S q (which becomes 
e l2 y/^ u ). In the phase estimation protocols that we pro¬ 
pose, the goal is to approximately measure the eigenval¬ 
ues of the unitary operator S p (and S q ). Once one knows 
the approximate eigenvalue of, say, S p , one could in prin¬ 
ciple apply an appropriate corrective displacement such 
that the eigenvalue of S p equals 1. In Appendix lAl we 
argue that it is not necessary to do so as any subspace 
characterized by fixed eigenvalues of S p and S q is a good 
code space (one can compare this Phase frame to a Pauli 
frame used in the description of stabilizer codes). 

Perfect code states which are eigenstates of S p and S q , 
can be parametrized 0 as basis vectors |zqv} (which 
form a complete orthonormal basis for the oscillator 
space) defined as 

\u,v) =e- lup e~ xvq \0), 

u £ (--A, Vtt], v £ VttT 2 ]- (7) 

Note that one has S p \u,v) = e l2 ^ v \u, v) and Z \u, v) = 
gZ^/iFu | U) v y Q ne can observe that a state | u, v) in g-space 
is simply a sum of peaks where the location of the delta- 
functions is shifted by u. while the wavefunction at each 
peak obtains a complex phase given by exp(— ivq). One 
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can similarly define approximate shifted code states as 

\u, v, A) = e ~ iup e~ ivq |0) approx 

OO n. 

°c Y e- 27rA2t2 e- {q - 2t ^ 2/i2A2) e- lvq \q + u)dq 

t =~OO ^ 

( 8 ) 

For such an approximate shifted code state the Gaussian 
envelope is identical, but each peak within the Gaussian 
envelope has obtained a phasefactor e~ lvq (we omit 
an overall phase factor) and the peaks are located at 
the shifted positions q = u mod 2 v / 7t. Of course, the 
width of the Gaussian envelope with standard deviation 
a = 1 /A should be larger than the maximal shift, i.e. 
A <C 1 /s/tt ~ 0.56, so that all shifts lead to approximate 
code states within the same overall envelope. It is the 
goal of the phase estimation protocols in this paper 
to produce such approximate shifted code states |it, v, A). 

It has been shown in 0 that it is possible to do quan¬ 
tum error correction using perfect code states which have 
undergone some (coherent or stochastic) distribution of 
shift errors of strength |w|, |u| < v^/ 6 , assuming other¬ 
wise perfect linear optical circuits. In Appendix El we 
review this way of doing quantum error correction. We 
will use the value, corresponding to a phase un¬ 

certainty 59 of the stabilizer checks of at most n/3 (but 
a phase uncertainty of 7 r /6 for the measurement of X 
or Z), as a rough figure of merit to represent how well 
one can prepare a code state using phase estimation. We 
will refer to this threshold as the \pK /6 shift error thresh¬ 
old. The ‘Steane quantum error correction’ analyzed by 
Glancy and Knill already assumes the existence of an 
approximately-prepared code state for which we devise 
a phase estimation method in this paper. Through ap¬ 
proximate phase estimation one can perform quantum 
error corection, hence the Glancy-Knill method may not 
be necessary to use at all. 

More formally, we define an effective shift error rate 
Ferror 6 of a state which makes reference to the \5k/6 
threshold as follows. Any density matrix p, i.e. an ap¬ 
proximately prepared code state, can be written in the 
basis | u,v) [ll[, as 


1 — f^/0 6 du Jf^ 6 / 6 dv |c(it,u)| 2 . For the approximate 
code state in Eq. O EH has shown that with n ~ 22 
photons, A = 0.15, the probability for shift errors 
beyond the shift error threshold is at most 1 %. 

If the shift errors are due to independent processes and 
the probability distribution p U v,uv = PuPv then one can 
use the effective shift error rate in p or q alone, i.e. 
pVW 6 _ 1 _ rvV 6 J and P ^/ 6 

1 error, q — a J — y/^K/Q (J ' LL ru <aiiu i error,p — a J — y/n/Q UjU r v 

to estimate the total error P e m( 6 . 

Note that this error rate is only a rough figure of merit 
as it makes reference to the shift error threshold of \J~Tt /6 
which has no intrinsic meaning if quantum error correc¬ 
tion is performed through other means. For an approx¬ 
imate code state one can also estimate the X or Z er¬ 
ror probability. The X error probability is the proba¬ 
bility that the state | 0 ) approx is incorrectly identified as 
encoding | 1 ) by means of perfect quantum error correc¬ 
tion. This is equal to the probability P e m(f q as shifts 
with |it| > y / ^'/2 lead to the state being incorrectly iden¬ 
tified. In [4| this error probability was approximately 
upperbounded by 2A exp(— 7 t/ 4 A 2 ) which leads to a 1% 
error probability for A = 0.5. Naturally, this error prob¬ 
ability is more optimistic than the threshold value, 
but since we expect the presence of additional (shift) er¬ 
rors during all protocols, using the threshold value seems 
a reasonable figure of merit. 

In practice, the code states are realized through an ap¬ 
proximate noisy implementation of the phase estimation 
measurement of S p and S q . Phase estimation by finite 
means, see Section El will typically output an estimate 
6 for 6 such that Prob(|0 — 0\ < 5) > 1 — e for some e 
and 5. This means that with probability at least 1 — e, 
one has \9 — 9\ = 56 < 5 phase uncertainty, which in turn 
corresponds to 5v = < <5/(2 v / 7t) shift uncertainty or 

error with probability at least 1 — e. Note that when one 
measures the eigenvalue of an operator such as Z = e l '^ q 
with phase uncertainty 59 < 5 (with high probability), 
then this corresponds to having a shift errror of strength 
at most 5u < 5/^/n. 


D. Physical Sources of Errors as Shift Errors 


/ V™ rVii 

du / du' 

-y/n j —y/n 

rV™/ 2 rV^/ 2 

/ dv dv’p uv y v ’ \u,v) (u',v'\. (9) 

J — y/n/2 J — y/n/2 

The effective shift error rate P e m ( 6 is defined as 

1 ‘y/ix/Q ry/n/6 

-P e ™ 6 = 1 - / du dv p uv ,uv, (10) 

so that for a pure state \if) = 
ff^du Sd^/d/ 2 dv c(u,v)\u,v), one has P e ™( 6 = 


In this section we discuss how levels of decoherence or 
inaccuracy of an arbitrary nature on the oscillator can be 
expanded into shift errors in p and q of small strength, at 
least when this noise is acting on states with a sufficiently 
low numbers of photons. We will consider this for natural 
physical sources of noise and inaccuracy such as dephas¬ 
ing, photon loss and quartic (self-Kerr) interactions of 
the bosonic mode. 

Following 0], any operator E acting on a single bosonic 
mode can be expanded as E = f_ d' y 0 ( 7 ) 11 ( 7 ) with 
complex coefficients 0 ( 7 ), i.e. an expansion into linear 
combinations of translations . 0 ( 7 ) with complex 7 in 
phase space. One can formally obtain the coefficients 0 ( 7 ) 
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as c(j) = ±Tr(ED(-y)) = ^ f da(a\ED(-'y)\a) = 
7 ) where C^( A) is the symmetrically-ordered 
characteristic function of the operator E (TH, HH, whose 
Fourier transform is the Wigner function of the operator 
E. Here da = dRe(a)dIm(a). Thus in principle one can 
find an error expansion for any operator by evaluating 
c(y) and considering whether its support is concentrated 
on small values of 7 . For example, the identity opera¬ 
tor / has c( 7 ) = S 2 ( 7 ), a delta-function in the real and 
imaginary part of 7 . 

At a more intuitive level, one can understand how¬ 
ever that the extent to which an expansion into small 
shifts is warranted, should depend on the number of 
photons in the state. Consider the action of an un¬ 
desired rotation exp(— iSa'a) (or photon decay opera¬ 
tor exp(— Sa'a)) in phase space. For simplicity, we ap¬ 
ply it to a coherent state |a): it is clear that such ro¬ 
tation corresponds to a larger state-dependent transla¬ 
tion for larger a, i.e. exp(— iSa^a) |a) = |aexp(— iS)) oc 
_D(a(exp(— iS) — 1)) |a) ~ D(—ia5) |a) for S <C 1. If 
exp(— iSa'a) = f dj c^D^) were only supported on 7 
with | 7| 2 < /(<$) with some function /(d) which is inde¬ 
pendent of the number of photons, then this contradicts 
the fact that there are states |a) with n = \a\ 2 large 
enough to be undergoing a large displacement |a<5| > | 7 | 

0 -. 

This implies that one should consider expansions (of 
physical noise operators in terms of shifts) which assume 
an upper bound on the number of photons in the oscil¬ 
lator space: let us call this upper bound on the number 
of photons n max . In that case it is clear from the above 
example that one can only make an expansion in small 
shifts, of strength at most y/n, of exp(— iSa'a) when at 

1 /o 

least one satisfies the inequality |d|nm ax < yfn. A sim¬ 
ilar argument can be given for a quartic interaction of 
form exp(— ie(a'a) 2 ). For such quartic interaction one 
can write, assuming e <C 1, 

(a| exp(ie(a^a) 2 )a exp(— ie(a 1r a) 2 ) |a) 

« aexp(—j|a| 2 e)) exp(—|a| 2 e 2 /2). (11) 

Thus ignoring fluctuations, one has an average phase 
space rotation (non-linear phase shift) in addition to co¬ 
herent state amplitude decay. As before, such average 
rotation and decay is only expressible in terms of small 
shifts of strength less than y/n when at least one obeys 
the inequality | a | | e | ' | e 1 t n x 

We note that these conditions are necessary, but not 
here proved to be sufficient as we have only applied these 
operators to coherent states |a) (while for example a su¬ 
perposition of Fock states |vac) + |n max ) undergoes a rate 
of change proportional to en^ ax by exp(— ie(a)a) 2 )). 

Since the approximate code states are superpositions of 
displaced squeezed vacua, we believe however that these 
criteria do give a good indication of whether a small-shift 
expansion is warranted and also capture the strength of 
shift errors. For the proposed protocol in Section mu 
we will show that the expected process of photon loss 


during the preparation of the code state has sufficiently 
small error rate so as to be deeply into this small-shift 
error regime. We will also discuss the (unwanted) cav¬ 
ity nonlinearities in the effective transmon qubit cavity 
Hamiltonian given the physical strength of parameters in 
a possible experiment. 

1. Shift Error Expansion 

Here we describe a general method to expand an error 
operator E in terms of shift errors. The error opera¬ 
tor E can arise either from open system dynamics or an 
unwanted unitary transformation. General open system 
dynamics of the bosonic mode can be modeled by a Lind- 
blad equation of the form 

P = -*[-ffideal(i) +V(t),p] + D(y/~Ka)(p) + V{yfya} a) (p ), 

( 12 ) 

with the compactly-defined superoperator V(X)(p) = 
XpX^ — ^(X^Xp + pX^X). Here -ffideai(i) is some ideal 
(time-dependent) dynamics (the execution of the phase 
estimation protocol of S p say, or some gate implemen¬ 
tation) and V ( t ) a possible perturbation or correction 
(e.g. an unwanted self-Kerr term K(a^a) 2 ). The Lind- 
blad equation can model photon loss D{yfna) of rate n 
or dephasing D(^a'a) of rate 7 . For simplicity, one can 
consider the effects of the different sources of noise sep¬ 
arately (and ignore Hideai(i) + V(t)) and construct a su¬ 
peroperator for the evolution in a short time interval, as 
S(p(t+r)) = JA Aip(t)A\. For the process of photon loss 
one has A 0 = I — a + 0((r«;) 2 n^ lax ) and Ai = y/nra. 

For dephasing one has Aq = I —^{a)a) 2 + 0({T"f) 2 n^ nax ) 
and A\ = yfyrafa. These Kraus operators are examples 
of error operators E. It is clear that one can express 
such operators as a low-order polynomials in a and 
and we will describe how such a low-order polynomial 
can be expanded in shift errors. 

Similarly, a unitary error operator E can be an over¬ 
rotation exp(— ida^a) or a nonlinearity exp(— ie(a'a) 2 ) 
where 6, e are given by an error strength times a time- 
scale. Assuming that Sn max C 1 or en 2 lax <C 1, one can 
Taylor expand this in terms of / and low powers of a and 
a). 

We now write a simple Taylor expansion in terms of 
correctable and uncorrectable shif t err ors for the oper ator 
\/~8a ~ Xi + ix 2 with x\ = p\J5/2 and X 2 = qy/5/2 
assuming 5n max 1. This expansion and its hermitian 
conjugate can then be used to expand operators such as 
a 2 ,a^a etc. With the Taylor expansion of arcsin(:r) ss 
x + + 0(x 9 ), we can expand x\ + 1 x 2 in 

terms of sin(xi),sin 3 (xi) and sin(x 2 ),sin 3 (x 2 ) etc. For 
example, in lowest-order one has 

VS a = sin [py/sj2) + sin 3 (py/S/ 2)/6 

+ i(sin(qy/6j2) + sin 3 (qy/5j2)/6) + 0 ((dn max ) 5/2 ), 

(13) 
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where the sin (p^/Sj2) and sin 3 (p^/Sj2) function (and 
similarly sm(qy/S/2)) can be expanded in terms of 

e ±i Ps /sj2 t e ±2i v ^72 j e ±3iy/s/2 eta In general we can 

thus write E = y/Sa = 7f correc tabie T A uncorrec t a bie 
where -^correctable contains all terms in the Taylor ex¬ 
pansion up to odd fc-th order such that kyj 5/2 ss 
\J 7t/2 and ||-K unc0 rrectable|| ^ O ((d/lnia.x) A )• We Can 
use this to expand an operator ( \/Sa) p (or (a)a) p / 2 or 

(id ) P ) aS (/\/~6(J.y > = -^correctable T -^uncorrectable where 

11 T^uncorrectable 11 — ^((^Imax )i +fe / 2 ) with k\f5 « \/tt/p. 

This shows that on an encoded state the amplitude of 
an uncorrectable error can be reduced by several orders 
of magnitude , depending on the strength of the errors. 

One thing to note is that in order to systematically ex¬ 
pand noise processes in terms of shift errors one should 
develop the full Kraus error operator or undesired rota¬ 
tion to a certain order in <5n max . In other words, depend¬ 
ing on which order in <5n max one chooses, one includes 
higher-order terms in the expansion of the error Kraus 
operators Ai or the unitary. 

One can interpret the protection that the code offers 
with the following example. Assume that photon loss 
from a cavity occurs for some time t at rate n such that 
P = Ktn max < 1. Without encoding, the error operator 
E = y/lda produces an uncorrectable error with proba¬ 
bility ~ P. For a cat-state encoding (see Sec II Bl) this 
process is correctable, but for the cat encoding, a two- 
photon loss error operator E = nta 2 is uncorrectable (as 
it is proportional to Z). Hence uncorrectable errors hap¬ 
pen with probability ~ P 2 . 

If we expand the 2-photon loss error operator in 
terms of correctable and incorrectable shift errors, 
the probability for the uncorrectable term is of order 
(ftTn max ) 1+fc / 2 = P 1+k / 2 where k ss \ yj 7rTl ™ ax 7 Thus for 
small P <C 1 the probability for an uncorrectable error 
resulting from 2 -photon loss can be much reduced if one 
uses a code that can correct shift errors as compared to 
a code which can only correct single-photon loss events. 


II. PHASE ESTIMATION 

The measurement of the eigenvalue e l6 of a unitary op¬ 
erator U, and the simultaneous projection of the input 
state onto the corresponding eigenstate U \ife) = e l6 \ife), 
is called phase estimation for U. In our case we have 
U = S p = D{\J 27r), say. Phase estimation can be exe¬ 
cuted by repeatedly running a circuit of a general form 
depicted in Fig.[3]for varying k and phases tp. Many vari¬ 
ants of phase estimation exist (see e.g. a recent analysis 
in [2l|) depending how or whether one varies fc, and/or 
whether one allows tp to depend on earlier qubit measure¬ 
ments, so using feedback, and how one infers the phase 
from the information obtained from the sequence of qubit 
measurement outcomes. 

We will first consider standard phase estimation as it 
is used in Shor’s factoring algorithm and argue that this 


method does not give the kind of approximate code state 
that we are looking for, at least not when we wish to use 
low numbers of photons. The standard phase estimation 
has the same performance as Kitaev’s phase estimation 
described in [22j ; both require realizations of U 2 for in¬ 
creasing k. 

The best form of approximate phase estimation opti¬ 
mizes the accuracy on the phase given a mean number of 
photons n in the approximately prepared code state as we 
always work with a bounded number of photons. In this 
respect it can be noted that for the standard phase esti¬ 
mation protocol described below, the phase uncertainty 
scales as 86 ~ 2~ AI ~ -4=, that is, shot-noise limited. 

v n 

Protocols which can achieve Heisenberg-limited scaling, 
that is 66 ~ i, have been considered and analyzed in 
detail in [23l - l25| (with an improved accuracy analysis in 
[26(|h One scheme that can reach the Heisenberg-limit is 
a modified form of the Kitaev’s phase estimation in which 
the circuit for each controlled- D(2 k y/2Tr) in Fig[3] is re¬ 
peated a number of times depending on k and M [2f| . 
One can thus expect that this Heisenberg-limited scheme 
performs optimally, but we do not consider it here as the 
experimental realization of the simpler schemes will al¬ 
ready be demanding. When one allows for a larger num¬ 
ber of ancilla qubit rounds, one can expect that switching 
to such scheme which uses U 2 or U 4 etc. is better. 

In sections III Bl and III Cl we analyze two simple phase 
estimation protocols, repeated or non-adaptive phase es¬ 
timation and phase estimation by feedback or adaptive 
phase estimation. These schemes only use U and thus no 
increasingly large displacements (microwave power). 

A. Standard Phase Estimation 

Any method for estimating the phase by sequentially 
acquiring bits of information is information-theoretically 
bounded: it is clearly optimal if each acquired bit gives 
us, with certainty, one additional bit of the binary expan¬ 
sion of = 0.0102 • • *, starting with the most significant 
bit 0i etc. This is close to what the standard phase es¬ 
timation protocol, as it is invoked in Shor’s algorithm, 
achieves [271 1 • H is known that this phase estimation pro¬ 
tocol can be realized sequentially, see e.g. [28l - [30[ . by 
coupling a sequence of M single qubits to the input state 
and applying controlled- U 2 gates (k = 0 ,...) with the 
ancilla qubit as control. This can be understood by re¬ 
placing the Fourier transform in the standard phase es¬ 
timation by the semi-classical Fourier transform whose 
quantum circuit is one in which the first qubit is mea¬ 
sured and its outcome is used to apply single qubit gates 
on the remaning M — 1 qubits. The next qubit is mea¬ 
sured and again its outcome determines single qubit gates 
on the remaining M — 2 qubits etc., see Fig. 5.1 in [27| . 
In this way the whole procedure can be implemented se¬ 
quentially, in rounds with one active control qubit per 
round. We depict one round of such a phase estimation 
protocol for U = S p = D(\Z2 tt) in Fig. [3] We can con- 








sider the measurement operator M x that is applied to the 
input state, upon getting an M -bit outcome x £ {0,1} M , 
with U = D(y/ 27t), including the compensating displace¬ 
ments in each circuit in Fig. [3] which center the code state 
around the vacuum state. One obtains M x |^i npu t), with 

i _ _ _o M i 

Mx = 2M E e ~ 2 ” iXt ' 2 D(V2^t)D(-V2^^—) 

te{ o,i} M 

2 m " 1 

« Y e~ 2nixt/2M D(V2nt). (14) 

t=—2 M ~ 1 

(Here relates to the fact that we simply approximate 
(2 m — l)/2 ~ 2 m_1 ). For phase estimation one can 
prove that, using M = M + log 2 (| + ancilla qubits, 
we obtain the best possible M-bit-estimate of the bi¬ 
nary expansion of with probability of success at least 
1 — e (when one chooses M = M, P su CC ess = yr)- Let 
9 £ (— 7 r, 7 r] represent this M-bit estimate. Let x roU nd be 
the bitstring outcome x rounded off to its M-rnost signif¬ 
icant bits. For a^ound/S^ < 1/2, we choose 9 = 27 r )w und , 

while for x r 0 und/2 Af > 1/2, one takes 9 = 2 ' xx ^gf' ld — 2i r. 

Knowing the first M bits of 9/{2tt) leads to an error 
59 < 27t 2 _a/ . This means that one can prepare an ap¬ 
proximate S p eigenstate with P success > yh (for which 
M = M + 2) such that the shift error is less than the 
\/rr/6 shift error threshold (correspond to 7 r /3 phase er¬ 
ror, see Section m if we take M > 4. 

There are a few ways to see that this method of measur¬ 
ing the eigenvalue of S p is not well suited for the approxi¬ 
mation that we are seeking. The protocol requires the im¬ 
plementation of controlled-.D(2 fe \/27r) for k = 0,... M— 1, 
Fig. [3] This means that the number of photons in the cav¬ 
ity mode grows exponentially with the number of ancilla 
qubits used. More precisely, one can calculate the num¬ 
ber of photons in an approximate code state produced 
by applying phase estimation of S p onto the squeezed 
vacuum state in Eq. ([3]), \if x ) oc M x |sq.vac). One has 

_ (sq.vac| M\afaM x Isq.vac) 

Tlx = -T- 

(sq.vac| MlM x jsq.vac) 

2-M-i 

E fc2 +^q 

fc= _ 2 M—i 

22 M —2 2 M_1 i \ 


where n sq is the expected number of photons of the 
squeezed vacuum state jsq.vac). Here we have used 
that D(—a)aD(a) = a + a and (sq.vac| D(y/2n(t — 
t')) jsq.vac) ~ 0 for \t — t'\ > 1 (the numerical value 
is 0(1CP 13 ) for A = 0.2). In principle n x depends 
on the outcome x , but the fluctuations with x van¬ 
ish when we approximate overlaps between different dis¬ 
placed squeezed states by 0. To leading order in M we 


have n x ~ 27r^— \-n sq ,, scaling exponentially with M 
as expected. For M = 4, one has n ~ 134 + n sq . This 
number is very high compared to the 22 photons in a 
Gaussian-approximate code state with A = 0.15 which 
has high probability to be below the shift error threshold. 

Another way to understand this is to consider the form 
of the measurement operator in Eq. m • Equating x ~ 

Ground, we have M x ^g oc J ] 2 = _ 2M _ 1 e~ lSt D(y/2nt). One 
can compare this measurement operator to the projector 
onto the space with S p = e ld which equals 

OO 

n Sp =exp( ie) = E e~ m D(V2^t), (15) 

t =—OO 

and note the similarity between the two operators. This 
shows that the filter used in standard phase estimation 
is a wide box-car filter with a hard cut-off, while we are 
seeking a smooth Gaussian filter. If we were to choose 
controlled-displacement operators with exponentially- 
growing displacements, one has to show that the duration 
of these gates is not necessarily exponentially-increasing 
as such growth in time would lead to more loss and de¬ 
coherence during the QEC cycle and code state prepa¬ 
ration. In principle large controlled-displacements can 
be obtained by large microwave power, see Appendix iDl 
but one may expect that the inaccuracy and undesired 
side-effects of exponentially-large displacements are also 
exponentially increasing. 

B. Phase Estimation by Repetition (Nonadaptive) 

The simplest way to estimate 9 £ (— 7 r, 7 r] of a unitary 
operator U \i/je) = e l6 \ipe) is to repeat the quantum cir¬ 
cuit in Fig. [4] with ip = 0 and (p = 7t/2. For the ip = 0 
measurement, each ancilla qubit then has a probability 
for outcome 0 equal to P v =o(0\9) = i(l + cos(0)) while 
for the p = 7 r /2 measurement one has P v= 7 r / 2 (O|0) = 
^(1 — sin(d)). Note that a simple repetition of the p = 0 
measurement is insufficient since Pm=o(O| 0 ) is the same 
for 9 and — 9. Thus one chooses p = 0 for half of the 
number of rounds/ancilla qubits and p = 7t/2 for the 
other half p 8 j . 

In an adaptive phase estimation, see Section III Cl one 
takes into account that the sensitivity of the probability 
distribution P v (0\9) = i(l + cos(d + p)) to 9 , i.e. 

dP *d6^ ’ a function of 9. One would like to optimize 
this sensitivity by choosing values for the feedback phase 
p which depend on previous measurement outcomes. In 
this optimization one chooses a next phase p such that 
the measurement with that p maximizes the sharpness 
of the resulting inferred probability distribution over 9. 
This method has been analyzed and described in detail 

in m, mi. 

One round of phase estimation for U = S p is depicted 
in Fig. [5j Instead of starting with an eigenstate, one 










9 


storage cavity 
qubit |0) 



FIG. 3. Phase estimation for the unitary operator S p = e~ l2 '^ p = D(y/2n) where D(a) is the displacement operator. Assume 
that phase estimation uses M ancillas prepared in |00 ... 0). The circuit in this Figure is repeatedly executed for k = M— 1,..., 0 
starting at k = M — 1. The phase ip in the single qubit rotation around the 2 -axis, diag(l, e* v ), will depend on the outcomes 
of all the previously measured ancillas. This sequential realization of phase estimation is identical to normal phase estimation 
as it merely uses a semi-classical realization of the quantum Fourier transform. Note that the circuit is identical to one in 
which diag(l,e 1¥> ) is moved before the controlled-displacement gate, which is the form of the quantum circuit in [2^]. Prior to 
the controlled-displacement gate, the cavity is (unconditionally) displaced so that the code states are symmetrically centered 
around the vacuum state and we minimize the total number of photons. 


IV’fl) 

qubit |0) 



li’e) 
qubit |0) 



FIG. 4. Phase estimation for a unitary operator U, with varying (feedback) phase <p, and U \ipe) = e lS | ipg). For repeated 
phase estimation one chooses either ip = 0 and n/2 for half the number of rounds. In a simple form of phase estimation with 
feedback the phase p is changed after each outcome so as to optimize the sensitivity of the probability distribution with respect 
to the currently estimated value for 8. 


would like to understand how one approximately projects 
onto an eigenstate of S p using this repeated measure¬ 
ment. We can first consider the action of this circuit 
on a coherent state |a) when the qubit is measured as 
x = 0,1: 

|a) -A |a - \Jn/2) + (-l) x e* v \a + \/it/2) , (16) 


has the simple form of a binomial sum of displacements: 
Mfc=o oc ^2 ( ^\D{V2n(m—M/2)) 

m—0 ' ' 

M/2 

« Y e~ 2t2/M D{V2^t). 

t=—M/2 


where the resulting state has not been normalized. Rep¬ 
etition of such a circuit on the resulting output state will 
thus produce a sum of coherent states on a line, each with 
a phase which depends on the feedback phase and the 
measurement outcome. If all phases add constructively, 
this distribution of amplitudes/weights of coherent states 
will be as in a Pascal triangle, hence binomial, as each co¬ 
herent state gets split in two equidistant coherent states 
at every round, see e.g. Fig.QJJa) in Appendix ITTl 

A more formal way of showing that the filter of this 
protocol is binomial and thus approximately Gaussian is 
as follows. Consider the measurement operator for M 
rounds with <p = 0. An outcome bitstring x € {0,1} M 
will correspond to a measurement operator which can 
be labeled by the Hamming weight k = wh(x) as the 
order of the outcomes of Is and Os is irrelevant. The 
measurement operator equals 


M k oc (/ + D(V2^)) M ~ k (I - D(V2n)) k D(—My^n/2) 

k M-k 


= E E 

pi—0 P2—o 


M-k 

P2 


{-l) Pl D(V^(pi T p 2 M/2)). 


When all measurement outcomes Xi = 0 (from which one 
would conclude that 0 = 0) the measurement operator 


One can roughly identify this measurement operator with 
the Gaussian-filtered projection operator onto approxi¬ 
mate code state in Eq. (U) with A 2 ft: (note that 

the measurement operator inevitably has a hard cut-off 
while the Gaussian-filter does not). This shows that for 
A = 0.2 one can choose M = 25/7T « 8 with a much lower 
number of photons than in regular phase estimation. 

Instead of postselecting on this outcome, we want to 
use the measurement data to estimate the value of 0, 
and include the data from the ip = n/2 measurement. 
We can analyze the efficiency of this method using the 
Chernoff bound; this argument is explicitly developed in 
fillf . Section III El gives numerical details on the phase 
uncertainty and the number of photons. Assume that we 
use M/2 ancilla qubits for ip = 0 and M/2 qubits for 
ip = n/2 such that P^c^O) is estimated as P v= o and 
Pyj= tt/ 2(0) estimated as P v = n / 2 - The estimate 0 is 
chosen as 

0 = arg(2P v=0 - 1 - *(2P ¥>=7r/2 - 1)), 

0 € (—7T,7r]. (17) 

Using a Chernoff bound in both cases, Prob(|P — P| > 
S) < 2e~ 2S M and some further bounding arguments 
[23 | . one can show that Prob(<50 > § \Jf{M)/M) < 
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4exp(—3/(M)/16) for any function f(M ) < M for 
M > 1. This implies that the probability to prepare 
a code state with phase uncertainty below 7r/3 is at least 
1 — 4exp(—3M/16). This argument shows that the num¬ 
ber of ancillas should be at least M > 8 and the prob¬ 
ability of failure will then rapidly vanish. We can note 
that the phase uncertainty 59 ~ -^= for M rounds, each 
of which adds 0(1) photons to the state, hence scaling 
in the expected (shot-noise) way. 


C. Phase Estimation with Feedback (Adaptive) 


The uncertainty of a probability distribution P(9) over 
phases 9 can be measured by the Holevo phase variance 
defined as V[P(9))\ = 5'[P(0)] _2 — 1 with the sharpness 
S[P(6)] e [0,1] given by 


S[P(0)] = |(e* e >| 



d,0e iB P(0) 


(18) 


For a 5-function distribution in 9 , the variance is 0, while 
for a flat distribution V[P(9)\ —> oo. One can show (see 
e.g. 0 ) that for small phase variance V[P{9)\ ss A 2 (9) 
where A 2 {9) = ((9 — (0)) 2 ) is the usual variance. 


Assume that we execute the circuit in Fig. [4] for M 
rounds so that one gets an M -bit estimate 9 of the phase 
9 and let a; [to] be the m-bit string of 0 or 1 measurement 
outcomes after to rounds. The adapative phases will be 
set to the values <pi,... p m , rn = 1,..., M and the ques¬ 
tion is how these will be chosen to maximize the sharp¬ 
ness and thus minimize the Holevo phase variance. We 
assume no priori knowledge on the phase 9 , i.e. the initial 
probability distribution over 9 , P(9), is assumed to be a 
flat distribution. This is reasonable as the squeezed (in q) 
vacuum state to which the S p measurement is applied is a 
superposition over eigenstates with rather uniformly dis¬ 
tributed phases 9 of S p . If we assume no prior knowledge 
about 9 , one may as well choose the first phase p i = 0 
which is what we do. 

It can then be argued, see the self-contained analysis 
in Appendix O that one should choose the next phases 
p 2 , ■ ■ ■, Pm as follows, depending on the previous mea¬ 
surement outcomes: 


p m = arg max 

<p 



d9e l9 P v (x[m\\9) , 


(19) 


where the probability P cp (x[m]\9) is the probability of 
obtaining measurement outcomes x\ ,..., x m given an 
eigenstate | ipe). This probability has a simple expression 
as the measurement results of each round are indepen¬ 
dent, i.e. 

m m / 9 _i_ . \ 

P V=Vm( x [ m \\ e ) = Y[ P Vi(Xi\9) = [Jcos 2 ( — —+Xi^ ) 

i— 1 2=1 ' ' 

( 20 ) 


For a (small) number of measurements, say, M = 
1,..., 20, one can simply solve this expression for the op¬ 
timal values for p %, P 3 ,... numerically given all previous 
possible measurement outcomes and store these optimal 
values as a look-up Table which is what we have done. 

Given that one has obtained a final measurement 
record x[M], how does one choose an estimate for the 
phase 97 In [31] it is argued that the estimated value 9 
should be chosen as 

/ 7T /»7T 

d9e ie P(9\x[M}) = arg / d9e ie P{x[M]\9). 

-7r J — TV 

( 21 ) 

Here we have used that P(9\x[M]) = 

P{x[M)\9)P{9)/P{x[M)). Note that P(9) is a flat 
distribution and P(x[M]) = 2~ AI which do not influence 
the arg function. The probability P{x[M)\9) 1 which 
implicitly depends on the feedback phases p \,..., p m 
was given in Eq. (17U1) . If all feedback phases are set to 
0, then 9 will be estimated as 0 or tt as the function of 
which we take the argument is real: this will however not 
occur since ^ 2^0 after the first bit has been generated. 
If we do not use any feedback, the value of this estimate 
9 coincides with the estimate in Eq. GH). 

In Fig. [5] we plot the probability for obtaining a phase 
error 59 = \9 — 9\ given a fixed number of rounds M 1 
averaged over random input phases 9 and runs through 
the protocol for both the non-adapative phase estima¬ 
tion protocol and the adaptive phase estimation protocol. 
One important difference between the two schemes is that 
in the non-adaptive scheme with p = 0 and p = 7t/ 2 one 
has a total of (-^ + l) 2 possible outcomes, as the order 
of the outcomes for a fixed value of p does not affect the 
measurement operator which is applied. This is clearly 
not true for adaptive phase estimation. Setting p\ = 0, 
we have two possible values for p 2 = p7^ ,x\ = 0,1. 
Then, there are four possible values for p$ lX2 etc. as 
the optimization in Eq. (1191) depends on all previous out¬ 
comes. The final estimate for 9 in Eq. (1711) will depend on 
all these phases tp* 1 '" 1 '- 1 and 9 can thus take on 2 M pos¬ 
sible values. These arguments suggest that the adaptive 
protocol can give a more accurate estimate of 9. Fig. [G] 
shows indeed how adapative phase estimation (APE) out¬ 
performs such simple phase estimation by repetition (PE) 
for a small number of rounds M. This difference will be¬ 
come more pronounced for larger M, see e.g. 0 , but 
the improvement is relatively small here. 

In order to understand better how well these phase 
estimation methods project the input state onto an ap¬ 
proximate eigenstate, we explicitly numerically generate 
the states that are created through these protocols in 
Section Hi El 


D. Preparing the Input State: Two Methods 

One has two alternatives for preparing a code state 
such as |0) (or |+}) approximately. As we have argued, 
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FIG. 5. One round of phase estimation for the unitary operator S p = e~ l2p ^ = D(y/2n) where D(a) is the displacement 
operator. The phase e l¥> of the -Rz-qubit rotation can be either fixed to ip = 0 and 7t/2 or adaptively changed per round 
depending on previous outcomes as in phase estimation by feedback, Section III Cl Prior to the controlled-displacement gate, 
the cavity is (unconditionally) displaced so that the code states are symmetrically centered around the vacuum state and we 
minimize the total number of photons. 


P(S6>e) 



FIG. 6. (Color online) Total probability P(89 > e) with 86 = 
\6 — 6\ versus t for M = 4, 8 (averaged over 6 £ (— n, 7r]). APE 
is the adaptive phase estimation protocol described in this 
section. PE is the non-adapative protocol where one sets the 
feedback phase </? = 0 for M/2 rounds and = 7t/2 for M/2 
rounds. For 86 < 7r/3, one is below the shift error threshold. 


|T(q)| 2 |T(q)| 2 



(a) (b) 


FIG. 7. Absolute values of two examples of wave functions 
of code states which can be generated by non-adaptive phase 
estimation for M = 4 and A = 0.2. After the eigenvalue 
measurement, the states are shifted onto an approximate code 
state using the phase estimate. The left state (a) is obtained 
if the first and the second measurement of both U and iU 
yield the same result, the right state (b) is obtained if both 
results differ. 


one can start the protocol in a sufficiently squeezed vac¬ 
uum state |sq.vac) and approximately measure S p (or 
squeeze in the other quadrature and measure S q to pre¬ 
pare |+)). Since the squeezed vacuum state is an ap¬ 
proximate eigenstate of Z one produces an approximate 
|0) state. How many dB of squeezing does one need in 
order to get a good eigenstate of Z ? Let the squeezed 
vacuum state be |sq.vac.) = J dq a q |g), Eq. ([3]), with 
A = e~ r and a q = -^e r l‘ 1 er < ? f ? r / 2 . The probabil¬ 
ity for shift errors above the shift error threshold is 


given by P e ™ 6 q < f^%/ 6 dq \a q \ 2 = Erf(e r v / 7r/6). 
One sees in Fig. [8] that the approach to 100% success 
probability happens around 8 dB of squeezing. Squeez¬ 
ing of itinerant microwave fields, instead of a confined 
cavity field, can be achieved using a Josephson para¬ 
metric amplifier (JPA). For example, in |32| one ob¬ 
tains A 2 = (A<? sq ) 2 /(Ag vac ) 2 a: 12% corresponding to 
A ss 0.35 and squeezing around 4.1 dB. More recent work 
has achieved squeezing of at least 10 dB |[33]. 

A more recent proposal was considered in [34j : in this 
paper it is analyzed how a squeezed microwave drive can 
be used to produce a squeezed vacuum state as the sta¬ 
tionary state of the cavity field (under dissipative photon 
loss dynamics). 


p_ 



FIG. 8. (a): Upper bound on the effective error probability 
Terror = Pe rrorfq for a squeezed vacuum state in the squeezed 
quadrature, (b): The average number of photons in the 
squeezed vacuum state depending on the gain G in dB. 

An alternative is to start the entire protocol in the 
vacuum state of the cavity mode and execute the phase 
estimation protocol both for Z as well as S p . The approx¬ 
imate measurement of Z is thus also an effective means of 
producing a squeezed state from the vacuum state in the 
cavity mode. The phase uncertainty 86 in the estimate 
of the eigenvalue of Z needs to be at most 7r/6 in order 
to be below the shift error threshold (instead of 7 t/3), 
twice as small as compared to the S p measurement. This 
suggests that if the phase of S p is measured in M rounds, 
then the phase of Z should be measured in 4 M rounds in 
order to give an overall similar shift error contribution. 
As Z = D(iy/n/2) as compared to S p = D(V2tt) the 
average number of photons added per Z-round is 1 /4 of 
a S'p-round, hence the total contribution to the number 
of photons would be about the same. These arguments 
show that it may be most advantagous to start with a 
squeezed vacuum state instead of selecting the code state 
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by phase estimation of Z as the ^-measurement would 
be rather lengthly (at least 32 rounds for M > 8) dur¬ 
ing which the code state is also decohering. At the same 
time, the Z-phase estimation protocol by itself may be 
an interesting novel way to prepare a squeezed state. 


E. Numerical Analysis of Two Phase Estimation 
Schemes for S v 


We apply the nonadaptive and adapative phase estima¬ 
tion scheme for S p with M rounds to a squeezed vacuum 
state with parameter A in order to create an approxi¬ 


mate code state |0) a 


The state obtained through 


the sequence of measurements is a superposition of dis¬ 
placed squeezed vacuum states, each with a phase which 
depends on the feedback phase and the measurement out¬ 
come x[M\: 


1 M 

|T(a;[M])} = ~^= ^ Cj{x[M])D{V2^{j-M/2)) |sq.vac). 


3=0 


with normalization N 


M 


3=0 M 


,(x[M])| 5 


( 22 ) 

(the 

approximation in the normalization comes from 
(sq.vac| D(y/2n(t — t')) |sq.vac) w 0 for \t — t'\ > 1). 
Here the coefficients Cj(x[M ]) = ^ TlfceSj 

where the subsets Sj are subsets of j = |5j | indices cho¬ 
sen from 1 ,...,M (without repetition). There are (^) 
such subsets Sj. Also, ipk is the feedback phase for round 
k and Xk is the outcome bit of round k. From the phase 
estimation one infers a value 9 £ (—7r, 7r], upon which one 
corrects |T(x[M])) -A |0(x[M])) approx = e ivq \ 4-(x[M]]) 

with v = Note that the additional displacement to 
center the code state around the vacuum does not affect 
the eigenvalue of S p as it commutes with S p . 

A few examples of the wave functions of states that one 
can obtain are shown in Fig. |TJ These examples show that 
one is not always guaranteed to get a good wavefunction 
and the effective probability of error can vary depending 
on the measurement outcomes. Because the relation be¬ 
tween the effective error rate and the measurement result 
is known, it is possible to use heralding or effectively use 
some post-selection. If we know (by simulation) that the 
obtained code state is bad, one would repeat the proto¬ 
col. 

In Fig. [15] at the end of the paper we plot the total 
probability P to prepare a code state with effective er¬ 
ror rate P e ™( 6 p for the (adaptive) phase estimation. As 
the effective error rate due to finite squeezing and the 
effective error rate due to the phase estimation scheme 
are unrelated (one corresponds to shifts in q , while the 

other corresponds to shifts in p), only the part P e Yror 6 p 
due to the phase estimation scheme is shown. The ef¬ 
fective error rate including both effects can be computed 

usinp- P^/ 6 — p^/6 I pyf^l 6 _ py/n /6 py/n/6 ™ 

UOAii & 1 error — 1 error,p i 1 error,q 1 error,p-* error,q* -L nc 


error rate P 


VtF/6 


error, q 


is shown in Fig. [8] 


Both protocols show a moderate probability to obtain 
code states with low error rates. What is striking is that 
the adaptive version of the protocol is much more reliable 
than the nonadaptive version and produces code states 
with an effective error rate below 1% in up to 94% of the 
cases. 

It is noteworthy that in most cases, both protocols re¬ 
quire some heralding/post-selection, as there is a finite 
probability to prepare a code state with a large effec¬ 
tive error rate. It should also be clear that using phase 
estimation in this way is much superior to simply post- 
selecting on obtaining the outcome x\ = X 2 = ■ ■ ■ = 
x M = 0 which occurs with an exponentially small prob¬ 
ability in M. If phase estimation is used for quantum 
error correction, then post-selection is not an option, as 
we do not wish to throw away the encoded state. 


1. Number of Photons 

We also consider the mean number of photons and fluc¬ 
tuations thereof for the adaptive and nonadaptive phase 
estimation of S p , applied to a squeezed vacuum state 
with parameter A. The expected number of photons in 
|4/(x[M])) (which is only slightly different from the num¬ 
ber of photons in |0(x[A/]))) approximately equals 

1 M 

n(x[M]) « sinh 2 (ln(A))+— ^ \(j-M/2)cj(x[M])V2n\ 2 , 

3=0 

(23) 

where we again used that (sq.vac| D(y/2Tr(t — 
t')) |sq.vac) ss 0 for \t — t'\ > 1 and the normalization N 
has been given above. In Fig. [8] we plot the contribution 
from squeezing on the mean photon number. In Fig. [9] 
we plot the total photon number after the measurement 
of S p . There, we see (n)e = f d,0P(x{M]\0) n(x[M]) 
with n{x[M]) as in Eq. (1231) as a function of M. We also 
plot the standard-deviation yj ((n — (rT) 2 )6>)e• Note that 
this is not the standard-deviation er(n) of the state itself 
which, as we discussed in Section [T] scales as n. As the 
contributions from squeezing and the phase estimation 
cannot be separated for finite squeezing, we show the 
photon numbers for A = 0.2. Note that both mean 
value and deviation depend on A. One can observe that 
n ss Mn/2 + n squeeze : this is what we expect as in each 
round we symmetrically displace the input state fu rther 
out to the left and right by an amount a = \J 7t/ 2, 
increasing the expected photon number thus by 7r/2.The 
total expected number of photons in the cavity for 
M = 8 and A = 0.2 (8.3 dB) is below 25 such that 
n ± <r(n) is below 50. 


F. Usage of Adaptive Phase Estimation in 
Quantum Error Correction 

Adaptive phase estimation of S p and S q can in princi¬ 
ple be used for quantum error correction. One cycle of 
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(a) (b) 

FIG. 9. The average number of photons {n)g after (a) adap¬ 
tive phase estimation of S p and (b) non-adaptive phase esti¬ 
mation of S p using M ancilla qubits (rounds) with parameter 
A = 0.2. The fluctuations yj ((n — (n) 2 )g)g around the av¬ 
erage number of photons reflect the different measurement 
outcomes in the protocol. It is noteworthy that these fluctua¬ 
tions are much smaller for the APE than for the PE protocol. 

quantum error correction will then consist of phase es¬ 
timation of S p followed by phase estimation of S q , each 
taking a certain number of M rounds. Each time one 
starts measuring S p one has prior information about its 
phase which should be used in choosing the next feedback 
phases. If the entire QEC protocol is noiseless, then the 
choice for the feedback phase of the first round of the S p 
measurement is simply estimated using the previous feed¬ 
back phases and outcomes of the previous measurement 
of S p , as determined by Eq. (ITiJl) . It is clear that for such 
a noiseless protocol, the number of photons thus keeps 
increasing round by round while the state gets closer and 
closer to a true eigenstate of S p (and S q ). In the im¬ 
plementation, see Section uni the number of photons is 
limited and at high number of photons the interactions 
that are invoked to execute the protocol are no longer 
accurate. In addition, photon loss is occuring continu¬ 
ously from the cavity during the protocol. This means 
that for quantum error correction the choice of feedback 
phases should ideally take into account an error model of 
the dynamics during the S q measurement and how this 
dynamics changes the current estimate phase for S p . We 
leave such more complete analysis to future work. 


III. PROPOSAL FOR REALIZATION IN 
DISPERSIVE QUBIT-CAVITY SET-UP 

We consider the following physical set-up. A super¬ 
conducting transmon qubit is capacitively coupled to a 
2D or 3D microwave high-Q storage cavity, as well as a 
low-Q read-out cavity. The storage cavity will be used to 
prepare a code state and the read-out cavity will be used 
to measure the state of the qubit. This is the same set-up 
as the cat state code experiments in e.g. 0] and we will 
use values for physical parameters which are similar as 
in this setup, see Table |T| We assume that a particular- 
cavity mode a with resonance frequency ui r couples most 
strongly to the transmon qubit and neglect the interac¬ 


tions of the qubit with other cavity modes, as well as the 
coupling to all higher-energy levels beyond the states |0) 
and |1). 

We assume that the interaction between qubit and 
storage cavity mode is approximately described by a sim- 

^bare 

pie Jaynes-Cummings Hamiltonian Hjc =- \—Z + 

oyada + g(<T_cd + cr+a). In the dispersive regime when 
| < 1 (A = j g ^hg detuning), one can make 

a perturbative expansion in and derive an effective 
Hamiltonian which equals 

H eS = (uj r - xZ)a)a - ^ q Z + O > ( 24 ) 

with dispersive shift \ ~ ^ + 0(g 4 / A 3 ). Here ui q = 
w bare _|_ ^ w h ere w bare - g ^hg bare q U bit frequency when 
the qubit is uncoupled to the cavity. The perturbative 
expansion is warranted when n < n cr it = A 2 / ( 4 g 2 ) where 
n is the mean number of photons in the cavity. The 
effective Hamiltonian shows that the resonant frequency 
of the cavity is shifted depending on the state |0) (+) 
or |1) (—) of the qubit, i.e. its frequency oj r —> oj p = 
uj r - FX- In this approximation we neglect a nonlinear term 
in the effective Hamiltonian of the form Z{a)a) 2 = 
x'Z(a^a) 2 , a nonlinear dispersive shift. 

In a more systematic approach the qubit-cavity Hamil¬ 
tonian can be obtained through first determining the nor¬ 
mal modes of the coupled LC circuits, after which the 
nonlinearity due to the Josephson junction is treated as 
a perturbation [ 35 | . This gives rise to an effective Hamil¬ 
tonian of the form 

iJ e ff = uir-a^ a+Cj q b^b—^^-(a) a) 2 -^^-(b^b) 2 — Xqr 0 ^ ab^b, 

( 25 ) 

x 2 

with Xrr = 2AT = -T- 32 —. Thus both qubit and cavity are 

4 Xqq 

represented as coupled nonlinear oscillators such that the 
anharmonicity for the qubit Xqq is relatively large. When 
one restricts the 6-oscillator to the lowest two levels, one 
can identify ui q « £o q — Xqq /2 and Xqr ~ 2%. Higher-order 
terms in Eq. (ESI) of the form {a)a) 2 tfb would describe 
the nonlinear dispersive shift. Note that in the model in 
Eq. the cavity self-Kerr nonlinearity is not present 
while its non-zero value has been determined in the ex¬ 
periments. 

The coupling of the qubit to the read-out cavity 
is described by a similar effective Hamiltonian as in 
Eq. CJ1) with a weaker dispersive coupling. 

In Table U we give the ranges of the relevant physical 
parameters and we choose several specific values in these 
ranges to demonstrate how well the protocol can be exe¬ 
cuted. We will take x/27r = 2.5 MHz and let the detuning 
be A/27r = 1 GHz in which case n cr it ~ ^ = 100. This 
upper limit on the number of photons is sufficient for cre¬ 
ating good code states: for an 8-round protocol M = 8, 
n « 25, see Fig.(9j so that n±cr(n) ;$ 50 (with a(n) ~ n), 
well below this limit. For this choice of parameters 
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g/2n ss 50 MHz. The strength of the neglected nonlinear 
term in Eq. (l24l) then equals /( 2n) « 6 ’^ a = 10 kHz. 
In the experiment described in 1161 ] in which cat states are 
created with about 55 photons, a value for xV(27 t) rs 4.2 
kHz is estimated for a detuning A/27T = 1.18 GHz and 
x/27t = 2.4 MHz which is considerably lower. The 
strength of the self-Kerr nonlinearity A'(ada) 2 in this ex¬ 
periment is estimated as K/(2n) « 3.61 kHz. The au¬ 
thors estimate n cr it < 300 which is a bit more optimistic 
than the estimate above. 

One round of phase estimation will last a total time 
Ground- During this time a quantum circuit depicted in 
Fig. [5] is to be executed. The protocol of a round consists 
of a short (O(lO)ns) interval T prep in which qubit (and 
storage cavity) are prepared. For the qubit this means 
it is put in the state ^=(|0) + |1}) (by a Hadamard of 
R x (tt/2) gate). After T prep , qubit and storage cavity 
should be coupled by a controlled-displacement trans¬ 
formation or a D{Z^f n/2) gate (in phase estimation for 
S p ). For this gate there are two options, as has been 
discussed in [l6j, [38] ■ 

In some physical set-ups the dispersive cavity-qubit 
coupling (both storage and read-out cavity) is not tun¬ 
able and is thus always ‘on’. Such a setup is non-ideal 
in various ways. When the dispersive coupling is al¬ 
ways on, it means that one should prepare the code 
states in a rotating frame (not the lab frame) which 
depends on the qubit state. For example, one chooses 
the frame in which the qubit is in the state |0) as the 
frame in which the cavity state should be unchanging, 
stationary (and we also look at the qubit in its rotat¬ 
ing frame as gates on the qubit are done relative to 
that). In this frame the effective Hamiltonian in Eq. CH) 
equals H e g = 2x|l) (l|a^a. It is important that during 
the qubit measurement, which takes up a considerable 
amount of time, this effective Hamiltonian induces no 
further rotations on the partially-prepared code state (or 
if it induces rotations, one should know what they are). 
This means that during qubit measurement the qubit has 
to quickly be reset to 10), in order to induce no further 
rotational dynamics on the cavity state. 

Another disadvantage of using a non-tunable \ is that 
the accuracy of single qubit rotations depends on the 
number of photons in the cavity. The qubit frequency 
given a storage cavity with n photons is, by Eq. m, 
given by ui q + 2xn. A microvave pulse which should 
rotate the qubit independent of the number of photons 
in the cavity should thus qualitatively take at least time 
T pu i S e = 2/W [ifj with frequency width W > 2ycr(n). 
This assumes that one sets the center frequency of the 
pulse at uj q + 2\n. Here <r(n) is the standard deviation 
in a code state, given in Eq. ©> which scales with 
n. In [l6} it was argued that the unwanted entangling 
of qubit and cavity due to single qubit rotations is a 
leading source of inaccuracies when one goes to higher 
photon numbers. Even though the number of photons 
in our proposed protocol for M = 8 will never be larger 


than 50, the protocol consists of many more single qubit 
gates than the experiments in [l(| so that these errors 
will accumulate. With a non-tunable x, the dispersive 
coupling to the read-out cavity is of course continuously 
on during the entire round while in the Figures it is 
suggested that measurement only occurs at the end of 
the protocol. It is understood that the coupling between 
transmon qubit and read-out cavity mode is smaller 
than the coupling between transmon qubit and storage 
cavity, for example in [l5j ] the dispersive coupling x/27r 
to the readout cavity was estimated as 0.930 MHz. At 
the same time, this coupling needs to sufficiently strong 
to provide a relatively short measurement time for the 
qubit as T roun( j will be largely determined by the length 
of the controlled-displacement transformation and the 
qubit measurement time. 

A third disadvantage of a non-tunable x is the the cavity 
Kerr nonlinearity which is present in Eq. (12511 due to the 
linear coupling between the LC oscillators: the ‘cavity 
mode’ is in fact a ‘dressed’ cavity mode which sees the 
Josephson nonlinearity. If x is turned to a small value, 
then this Kerr nonlinearity will be correspondingly small. 

For these reasons, we imagine that the dispersive cou¬ 
pling x is tunable and can be turned ‘on and off’ during 
the execution of a round. Thus during each round, the 
dispersive storage cavity-qubit coupling is ‘off’ (i.e. 
to a low value) during qubit preparation and single qubit 
rotation (lasting time T prep ) and possibly unconditional 
cavity displacements (of the form D{—\Jtt/2)). Then 
X is rapidly turned on after T prep and stays on for a 
time-interval T during which the controlled-displacement 
gate in Fig. [5] acts. The coupling x is again rapidly 
turned off during the last interval T post during which 
the qubit state undergoes further single qubit rotations 
and is being measured. One can similarly imagine that 
only during the measurement time the coupling to the 
read-out cavity is turned on. The turning off and on 
of x could be achieved in two possible ways. One can 
increase the cavity-qubit detuning A (and hence reduce 
X') by altering the qubit resonant frequency w^ are using 
a flux-tunable transmon qubit; such switching can take 
place in 0(1) nsec. An alternative is to have a variable 
qubit-cavity capacitive coupling {4lJ, |42[ which is turned 
on and off in a 0(1) nsec time window. 

Let us now discuss the realization of the controlled- 
displacement gate. Note that if x is in principle turned 
‘off’ unless the cavity state needs to be manipulated, 
it means that we are preparing the code states in the 
rotating frame of the cavity at frequency uj r (while 
single qubit gates are performed in the frame rotat¬ 
ing at w£ are ). When the coupling is turned on, the 
effective Hamiltonian in these rotating frames is then 
H e s = —\Za^a — \\Z (neglecting nonlinearities). 
The additional Z-rotation in this Hamiltonian leads 
to a phase accumulation on the prepared qubit state 
^(|0) + |1)) which should be taken into account when 
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Transmon qubit and bare cavity frequency ^3 — 11 GHz 


Qubit T 1 /T 2 time 

3D (storage) cavity lifetime T cav 

2D CPW (storage) cavity lifetime T cav 

Controlled-displacement pulse time T pv ,i se 

Dispersive shift x/2n 

Qubit measurement time tmeas 

Single qubit gate 


10 — 100 /rsec 

55/rsec [15], 1 msec (Q > 10 7 ) [36] 
200 /rsec (Q > 10 6 ) [37] 

25 - 100 nsec [38J 

1 - 20 MHz 

200 - 300 nsec [15, 39] 

5—10 nsec 


TABLE I. Ranges of some relevant parameters. The quality factor values (and lifetimes) represent the internal losses inside 
the cavity while in our protocol the total quality factor of the cavity Qtot, including its intended but flexible coupling to the 
outside world, is the relevant parameter. 


considering what single qubit Z-rotation of the form 
diag(l, exp(f<^)) is done during T post . 

In principle, one can enact a controlled-displacement 
gate by supplementing this dynamics by driving the cav¬ 
ity with a microwave pulse. During the time interval T, 
a microwave drive of duration T pu i se « T is applied. The 
drive tone of this pulse uid = + y is resonant with 

the cavity mode when the qubit is in the state |1) but 
off-resonant when the qubit is in state |0). Hence, one 
expects a large cavity displacement in the resonant case 
and a small neglible displacement when the drive-tone 
is off-resonant, thus enacting a controlled-displacement 
gate. Again, at an intuitive level, the frequency width 
W of the pulse should be sufficiently narrow as compared 
to 2\, which is the difference in resonance frequencies of 
the oscillator given the qubit state, so that the pulse has 
few photons at frequency u> r — x- If we center the pulse 
at u> r + x, then for a Gaussian pulse with T pu i se = 2/W, 
one should at least require that 2% > ^ = 1/T pu i se or 
X7 pu i se >1/2. This heuristic argument assumes that the 
cavity decay rate k which sets the width of the resonance 
obeys k <C W which is warranted for high-Q storage cav¬ 
ities. 

During this pulse, the qubit-cavity coupling — x^ a ^ a 
also induces a qubit-state dependent phase space rotation 
on the cavity state, i.e. of the form exp (ixZa) at) : which 
is in fact undesired , see Fig. [T3]in Appendix iDl It is thus 
important to choose the interaction-time T such that this 
difference in rotation angle is 27 t which implies that the 
newly created superposition of displaced squeezed states 
stays on a line. This warrants a precise choice for T 
namely x'T = ir. In the analysis in Appendix[D]we model 
the drive and the corresponding displacements and ro¬ 
tations assuming a simple square pulse (with zero rise 
time). The conclusion of this analysis is much less nega¬ 
tive than what has been stated in [38| where it is written 
that a direct controlled-displacement gate requires a time 
T pu i S e ~ 30/x which would be quite long. 

The alternative is to do a controlled-displacement gate 
through a sequence of two controlled-rotations inter¬ 
spersed with an unconditional displacement [43], explic¬ 
itly shown in Fig. [TU] 

Let R(9 ) = exp(— iOa^a) be a rotation in phase space. 
The conditional rotation R(—Ztt/ 2) = exp^waZn/2) 


is obtained by TT e ff, he. coupling the transmon qubit 
for time t with tx = 7r/2 to the cavity mode. The 
inverse rotation equals R(Z-k/ 2) = XR(—Zir/2)X 

where X is a Pauli X (7r-rotation) on the qubit. The 
advantage of this way of doing a controlled-displacement 
is that its implementation is always fast, as it requires 
0(10) nsec. unconditional displacements and single 
qubit-rotations as well as a total qubit-cavity interaction 
time of T = n/x during which the cavity state is rotated 
depending on the qubit state. For our chosen x, one has 
T = 7r/x = 200 nsec. Note that the last X-gate on the 
qubit can be absorbed in the post-processing stage dur¬ 
ing which one thus applies a single qubit rotation which 
is composed of Pauli X in Fig. [TU] the possible feedback 
rotation diag(l, exp(?'<^)), including corrections for 
accumulated phases, and then a Hadamard gate, shown 
in Fig. [5] However, in order for this implementation 
to work, one needs to turn off the dispersive coupling 
X in the middle period during which an unconditional 
displacement and a single qubit rotation act. Thus this 
scheme would require one to turn x on (and off) twice 
during a round, while a direct conditional displacement 
would allow one to turn x on and off once. 

In either implementation of the controlled- 
displacement gate, T roun( j will be largely determined 
by the duration of the controlled-displacement gate 
and the qubit measurement time. If we assume 
that the qubit measurement time is approximately 
300 nsec, we have T roU nd ~ 500 nsec. During this 
time qubit decoherence (see Table [D is negligible as 
T roU nd/Ti(= 50^sec) = 1 x 10 -2 . If we assume a 3D 
storage cavity lifetime of 1msec (see Table [J), then 
during the execution of M = 8 rounds for the full phase 
estimation of S p , in total lasting Tpg = 4/zsec, one 
has TpE/T cav = 4 x 10~ 3 . One cycle of quantum error 
correction using the measurement of S p and S q will thus 
last Tqec— cycle = 8 /rsec, still considerably faster than 
the cavity decay time. 

For 3D cavities, no flux-tunable transmon qubits have 
been used so far. For a 2D co-planar waveguide cavity, 
to which transmon qubits have been coupled with a 
tunable x> the best lifetime can be 200/zsec (see Table) 
which would give estimates which are factor of 1/5 worse 
as compared to a 3D cavity. 


cavity 

qubit 
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FIG. 10. Realization of a controlled-displacement gate via two controlled-rotations for one round of phase estimation of S p 
with q = \/27r or S q with a = —i\J 27r, see Fig. [5] Here R(—Zn/2) = exp(icdaZ 7 r/ 2 ). 


The adaptive phase estimation protocol uses feedback 
which implies that the phase in a single qubit rotation 
in the next 500 nsec round depends on the outcome of 
the qubit measurement in the previous round. In current 
technology it is possible for a qubit measurement out¬ 
come to determine the execution of a single qubit gate 
200 nsec later, so this does not pose a problem in the 
implementation. 

A. Shift Error and Noise During The Protocol 

There are two sources of errors during a round of the 
protocol, namely errors on the transmon qubit and direct 
errors on the cavity mode. We first consider errors on the 
transmon qubit which can propagate to the cavity mode 
in two ways. 

We can call a protocol strictly fault-tolerant if an error 
with low probability (or amplitude) on the qubit can in¬ 
duce a shift error u or v with low probability/amplitude 
and low strength |w|, |u| -C y ^/2 on the cavity mode. 

An error on the qubit during a round can alter the 
measurement outcome of this qubit (for example an error 
which flips the measurement outcomes). This will lead to 
an incorrect estimate for 9. However, if the rate at which 
these errors take place is sufficiently low, so most rounds 
are error-free, then it should be clear that the error on 
the estimate is also small, hence the induced shift error is 
small, implying fault-tolerance. One can see this as fault- 
tolerance which arises from the repetition of the rounds; 
how well this works will depend on whether the rate of 
qubit errors is sufficiently low. 

Another type of error on the qubit can directly result 
in a shift error on the cavity mode. For example, dur¬ 
ing the execution of the direct controlled-displacement 
or rotation, the qubit state decays to | 0 ), resulting in 
a displacement error which can be as large as D(y/2n). 
Such error shows that the phase estimation protocol is 
not strictly fault-tolerant, since, even though this error 
is very unlikely, it does lead to a large incorrectable shift 
error. Using several, say k , ancilla qubits per round 
and putting them into a cat state |00 ... 0 ) + 111 ... 1 ) 
so that each one of the k qubits is only used to perform 
a much smaller controlled displacement D(y/2n/k) can 
mitigate this problem and make the protocol more fault- 
tolerant. Effectively, we are trading a low probability for 
a large shift error to a larger probability (as one uses 
more qubits) to have only small correctable shift errors. 

Given the fact that the accuracy of single qubit 


gates and controlled-rotations is quite high and qubit 
decoherence on the time-scale of a round very low, we do 
not anticipate that qubit errors are the dominant source 
of errors. 

Let us next consider direct errors on the cavity mode 
which can result from photon loss, a self-Kerr nonlin¬ 
earity or a nonlinear dispersive shift. In principle, for 
short enough time-intervals all these processes can be ex¬ 
panded as linear combinations of shift errors and these 
shift errors can be propagated through the ideal circuit 
in a round, or in between rounds, such that they remain 
shift errors which gradually add up in strength. This is 
true as the ideal circuit only contains (conditional) dis¬ 
placements and (conditional) rotations. 

Given the points made in Section IID 11 one can con¬ 
sider the effect of cavity decay during preparation and 
the presence of nonlinearities which we expect to be the 
dominant source of shift errors. For cavity decay with 
rate n with k = 1/T cav of a 3D cavity, one can easily 
meet the condition 6 = kTpe ~ 4 x 10 ~ 3 -C —— when 

f^max 

Umax created during the 8 -round protocol is no more than 
50 photons. The arguments in Section IIP 11 then im¬ 
ply that cavity decay modeled by some Lindblad equa¬ 
tion induces dimensionless shift errors u, v of strength 
| it |, |u| < 0 (\/kT pe) ps 0.06 (using only the first term in 
the shift expansion in Eq. m during the phase estima¬ 
tion measurement of S p . 

We can also consider the effect of the self-Kerr non¬ 
linearity in Eq. (12511 or the nonlinear dispersive shift 
x'Z(a^a) 2 where the strength of x' /(27t) and K/(2tt) can 
be taken to be, say, 4 kHz. If we represent the overall 
nonlinearity during the protocol as exp(—*e(cda) 2 ) one 
has |e| = KTpe = 27r x 16 x 10 -3 ps 0.1 which is two or¬ 
ders of magnitude larger than — r — for n max = 25. This 
means that the shift errors due to the nonlinearity will 
add up throughout the protocol and do not necessarily 
remain small correctable shift errors. 

Even though the physical values of these numbers are 
not clad in stone, there is a simple argument which 
demonstrates that higher-order terms in a pertubatively- 
derived dispersive coupling Hamiltonian will lead to er¬ 
rors which require a more thorough analysis. This argu¬ 
ment comes about from the fact that the running time of 
the protocol, let’s say T roun( j, is determined by x ~ 5 2 /A, 

i.e. T r ound ~ 1/y. The next order in the perturbative 

4 

coupling between cavity and qubit scales as ^3 and we 
want that its overall effect during T rounc j to be sufficiently 
small, but still have a sizeable number of photons in the 
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cavity. We consider 


e \n„ 


^3 7roundn max 1 


a 2 2 

- 71 

^2 '‘'max 


n„ 




< 1, 


(26) 

while rimax should be 10 — 100 photons in order to ob¬ 
tain a good code state. One has n max < n cr it but the 
stringent condition in Eq. m would require making ?z cr it 
much larger. Since n cr it should be enlarged without mak¬ 
ing the protocol last much longer (as this would enhance 
the strength of other sources of errors), it would be best 
to work at a larger detuning A rather than a smaller 
capacitive coupling between transmon qubit and cavity 
(coupling g). 

One should note that this problem rapidily gets better 
at higher-order terms in perturbation theory. Assume 
that the first unwanted term comes in at fcth^ order in 
perturbation theory, i.e. with strength g (^) , so that 

Eround*? (a)^ 1 n max <C 1. One then obtains the condi¬ 
tion 


4n cr i t 


f ^max \ 
\ 4^-crit / 


k/2 


< 1. 


(27) 


For a quartic term k = 4 this condition is quite de¬ 
manding but for the next higher order term k = 6 and 
n m ax <C 4n cr it the condition is much more mild. 

These arguments suggest that it is better to not treat the 
nonlinear dynamics as a source of errors, but rather treat 
it as a source of known systematic errors. One approach 
is to try to actively cancel their effect during the evolu¬ 
tion so as to obtain the same code states j44|. Another 
approach is to take as a given that these interactions ex¬ 
ist and seek a code formalism that captures their effect. 
Assuming that % is tunable, one may at least hope to 
reduce the nonlinear dispersive shift and the self-Kerr 
nonlinearity in strength during the time that the disper¬ 
sive interaction is off. 

These arguments show that further numerical analy¬ 
sis including the nonlinearity and open system dynamics 
modeled by a Lindblad equation are warranted to assess 
their overall effect during the phase estimation protocol. 


IV. DISCUSSION 


(sequentially) coupling a qubit to the oscillator to gather 
phase information, we are getting little information, at 
most one bit, per qubit used. At the same time, this 
transmon qubit is effectively realized as a non-linear os¬ 
cillator as in Eq. (1251) . So can one not get information 
faster by coupling the storage cavity to another oscilla¬ 
tor? 

Assume that we can prepare one ancilla oscillator in 
the state \q a = 0) where q a is its position quadrature 
(or some squeezed version thereof) and the other ‘stor¬ 
age’ oscillator is in some initial state \i/j s ) = f dp a p \p). 
Assume an interaction between the oscillators of the 
form Hi n t = Ps <8> Pa for some time t = 2y / 7r so that 
e -iH irit t (g) \q a = 0) = / dpa p | p) | q a = 2 py/n). If 
we can measure the q quadrature of the ancilla oscil¬ 
lator modulo 2-7 t, i.e. we determine gwas = qa mod 27r, 
we project the system oscillator into a superposition of 
states with p = ^=r mod spS. Shifting back this pro¬ 
jected state depending on g mea s gives the output state 
YlkL -oo a p\P = thus P re P ar ing an eigenstate of 

S p . A similar procedure for S q on this state would then 
prepare the full code state. The linear coupling between 
the two oscillators is very natural; superconducting LC 
circuits with capacitive coupling realize such interaction. 
However, measuring a quadrature modulo 2tt is not sim¬ 
ple: one has to make sure not to get too much informa¬ 
tion. If the ancilla oscillator were modified to include 
a Josephson junction, it would make one of the quadra¬ 
tures, namely the phase variable, 27r-periodic and thus 
suitable. A rapid turning on of the Josephson junction, 
and thus a rapid change of the quadratic potential in q or 
(j) to a periodic cosine potential U((j>) = —Ej cos(< j>) could 
freeze this state (taking Ej Ec , potential energy much 
larger than kinetic energy) and information could be read 
out. It is an open question whether a rapid turning-on 
of a Josephson junction (see e.g. M for a similar cir¬ 
cuit switch idea) is experimentally feasible and what the 
details of such a scheme would look like. 


V. ACKNOWLEDGEMENTS 

We would like to thank David DiVincenzo and Robert 
Schoelkopf for comments and feedback. We acknowledge 
support through the EU via the programme SCALEQIT. 


In order to implement phase estimation in a shorter 
amount of time, it is possible to couple several transmon 
qubits simultaneously to a 2D (or possibly 3D) cavity. It 
may be interesting to consider whether a form of Shor 
quantum error correction (using cat state ancillas) in¬ 
stead of Steane quantum error correction is possible for 
the GKP code states. It will be worthwhile to numeri¬ 
cally study the code state preparation protocol including 
the effect of nonlinearities and open system dynamics. 

A somewhat different scheme for preparing code states, 
eigenstates of S p and S q , involves two oscillators as fol¬ 
lows, borrowing an idea in pihj . The basic idea is that by 


Appendix A: Use of Phase or Displacement Frame 
and Logical Gates 

We prepare a code state or perform error correc¬ 
tion through phase estimation of the displacement op¬ 
erators S p and S q . The phase estimation proto¬ 
col outputs an estimate for the eigenvalue e l6p (e* 6 * 9 ) 
of S v (S q ), while simultaneously projecting the input 
state onto an approximate eigenstate \4’e p ,e q ) of S p 
and S q . Error correction would then correspond to 
displacing this oscillator state by the corrective dis- 
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placement Dg p g^ = e iqep ^ 2 ^ > e ipSq ^ 2 ^^ such that 
A? p ,e q \ipe p ,8 q ) ~ 1 ^g P =o,e q =o), an approximate code 
state. 

Here we show that it is not necessary to do these 
additional displacement operations, but one can work 
with a phase frame similar as the Pauli frame for 
stabilizer codes mi,!. We will use the notation A' = X 
and U for operators which have the action of Pauli X 
and a unitary gate U on the states in the code space. 
We call such operators logical operators or logical gates. 
In addition we write Xg p for the Pauli X operator on a 
code space labeled by a phase 9 P for S p . 

Let us first assume that we work with perfect phase 
estimation and perfect code states. The operator e~ l ' z * v 
has been called X because it maps the perfect state 
|0) onto |1) and vice versa. We can note that X is not 
hermitian, we have in fact A' = A S p . The action of A' 
and X ^ on the perfect code space is however identical, 
thus on the code space A' = A'^ so that X 2 = I. If we 
use a code space labeled by the eigenvalue e l6p of S p , we 

have A = e lBp I. This means that we need to redefine 
the logical X operator on the code space characterized 
by eigenvalue e l6p as Xg p = D\ g q XDg P} g q = Xe~ lBp ! 2 

-2 

for which X g = I on the 9 P code space. Similarly, 
z 8 q = D g p ,g q ZD 8 P ,e q - Clearly, Xg p and Zg q only differ 
from X and Z by phases so transform similarly. 

It is known what operations are necessary to perform 
logical gates on the code states in the 9 P = 0, 9 q = 0 code 
space, such as the Clifford gates (the CNOT gate, the 
Hadamard gate, the S gate) and the T gate, see !, as 
one can verify their proper action on the logical operators 
X and Z. 

The logical Clifford group gates can be realized by lin¬ 
ear optical transformations, i.e. linear transformations 
on the set of positions and momenta of n oscillators, 
(qiiPi, ■ ■ ■, qn,Pn), which preserve their commutation re¬ 
lations. A circuit for the CNOT gate is shown in Fig. [T2l 
The Hadamard gate represents a 7r/2 phase delay enact¬ 
ing q — > p and p -A —q. The S = diag(l,i) gate enacts 
the transformation q —> q and p -A p — q. 

Under the action of these gates, shift errors (dis¬ 
placements) remain shift errors, i.e. a general dis¬ 
placement operator on these n oscillators of the form 
D = exp (iJ2i=i( a iPi + PiQi)) (with real a u ft) trans¬ 
forms to D' = UDW where D' has coefficients WuP'i}- 
One can propagate the shift errors through a Clifford 
circuit: small shift errors in several modes can add 
up to large shifts in one mode, similar as Pauli errors 
can propagate to become high-weight incorrectable 
Pauli errors. Furthermore, even though the Clifford 
gates do not amplify shift errors, shift errors inside, 
say, the realization of a CNOT gate can get somewhat 
(de-) amplified in strength as the circuit uses squeezing 
and squeezing acts e.g. as p —> e r p 1 q —> e~ r q. 


For logical Clifford gates U one can verify their action 
on the 9 p ,9 q code space. One has 

U Wo p ,g q ) = UDl pe u'u |^o, o) = De p ,e q U |V>o,o), (Al) 

where Dg p: g p is the new shift correction (displacement). 
This means that one never needs to do the corrective 
displacement, but can just keep track of the phase frame 
in software, just as for the Pauli frame in a computation 
with only Clifford gates. 

When the encoding is only approximate (and we know 
the eigenvalues of S p and S q approximately), every 
code state is in principle given by a perfect code state 
and a distribution of shift errors, i.e. a density matrix 
expanded in the shift error basis as in Eq. ([9]). The 
gates are still defined by their action on the perfect 
code states. The distribution of shift errors or the 
approximate encoding may be slightly different for every 
code state when these are generated, say, by the protocol 
in this paper. It is important to note that this does 
affect the proper functioning of the Clifford group gates 
as these gates only propagate the shift errors between 
logical qubits. 

For quantum universality, one needs a T = 
diag(l, exp(«7r/4)) gate. This gate can in principle be 
realized using an ancilla T|+) or alternatively the an- 
cilla \H = 1) oc SHT |+) which is a +1 eigenstate of 
the Hadamard gate, in addition to Clifford group gates 
(see e.g. !). The logical Hadamard gate H: X o Z 
equals exp(f : |a^a) which means that a +1 eigenstate of 
this operator has a photon number which is 0 mod 4. 
One can create the state \H = 1) by doing phase estima¬ 
tion for S p and S q on the vacuum state |vac) for which 
e *f a+a | vac ) = |vac) and post-selecting on the outcomes 
S p « 1 and S q « 1. The resulting state is encoded and a 
+1 eigenstate of H as H commutes with S p and S q . One 
can imagine that noisy ancillas \H = 1) are further dis¬ 
tilled into fewer higher-quality ancillas before being used. 
It is important to note that a +1 eigenstate of H for the 
phase frame 9 q = 9 p = 0 is not related by a displacement 
to a +1 eigenstate of H in an arbitrary phase frame 9 qi 9 p . 
So before one uses such ancillas one has to at least know 
or physically fix the phase frame that is in use. We leave 
a realistic method to implement T (and Clifford) gates 
on the approximate code states in circuit-QED hardware 
to future work. 


Appendix B: Quantum Error Correction Using 
Encoded Ancillas 

In this section we briefly review a form of quantum 
error correction suggested in ! and a simplified imple¬ 
mentation in [lH . which, for stabilizer codes, is called 
Steane Error Correction. In Steane Error Correction the 
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errors are corrected by coupling the encoded data to an 
encoded ancilla using a CNOT gate and measuring the 
encoded ancilla. The advantage of using this procedure 
for quantum error correction is that the quantum cir¬ 
cuit is fault-tolerant, made from linear optical elements, 
phase shifters and squeezers and is fully deterministic. 
Fault-tolerance means that shift errors are not amplified 
in strength by the circuit (but they do propagate) as it 
involves only linear optical components. A disadvantage 
of the method is that it requires an ancilla code state 
which has to be first prepared with sufficient accuracy. 


e -iup 


|#) 

1 +) 


-•- 

— Z?( — q mod y/n)) 

■—© 

—-d) 


cilia state are of strength at most |w|, |v| < v^/6, errors 
remain correctable through the repeated application of 
these circuits, assuming that the linear optical circuits 
themselves are faultless. If the linear optical circuits are 
faulty, one has to propagate the shift errors that occur 
during the circuit forward which will result in a some¬ 
what lower threshold value. 


In 


Out 



-31.7° -31.7° 


FIG. 11. Steane Error Correction for shifts in q generated 
by the error shift operator e~ lup which are detected by the 
stabilizer check S q . The measurement of the ancilla in the 
encoded |+) state is a homodyne measurement of q. The 
displacement D(—q mod y/i r) is a correction which is like 
a Pauli frame for a stabilizer code and does not physically 
need to be realized. Note that the modulo function should 
be taken in the interval (—y^r/2, y/Tt/2]. In the main text we 
have referred to |+) as |+). 


The basic circuit for correcting shifts in q (detected 
by S q ) is shown in Fig. 1111 A similar circuit holds 
for S p where the encoded ancilla is prepared in |0) 
(denoted as |0) in the main text). The CNOT gate 
(a linear optical gate) has the following effect on the 
two quadratures of control (c) and target (t) mode: 


qc -> q c ,qt -> q c + qt,Pc -» Pc - Pt,Pt Pt- The CNOT 
gate cannot be implemented with only beamsplitters 
and phase-shifters as the symplectic transformation 
matrix of the CNOT, i.e. the matrix applied to the 
vector (q c ,Pc, qt,Pt) is not orthogonal, while it is always 
orthogonal for any passive linear optics transformation, 
see e.g. f47j. The CNOT gate can be realized by beam¬ 
splitters and squeezing [48}. The optimal CNOT circuit 
expressed in terms of these elementary components is 
depicted in Fig. [12] 


In the Glancy and Knill method in fll[ the CNOT 
gate which uses two beam-splitters and two squeezers 
is replaced by a single beam-splitter and a single (but 
stronger) squeezer. The y/Tt/6 shift error bound is arrived 
at by arguing how shift input errors on data and ancilla 
propagate through the quantum error correction circuit 
leading to shift errors which have to be corrected in the 
next round. 

Basically a shift error exp(— iup) exp(— ivq) that needs 
to be properly corrected in steady state (which requires 
|u|,|u| < V^/2) results from the propagation of three 
separate shift errors u = ui + u p + u q , v = vi + v p + v q , 
which are the initial errors w/, Vi and the errors that oc¬ 
cur somewhere in the quantum error correction circuits 
for S p (u p ,Vp) and S q {u q ,v q ), hence all |'Uj|, |i>j| < y/n/6. 
Thus if all initial shift errors on encoded data and an- 


FIG. 12. (Color online) A CNOT gate on two bosonic modes, 
target and control. S denotes squeezing and its strength in dB 
(see squeezing convention below Eq. (|4|l). The phase shifters 
shown explicitly can of course be absorbed in the action of 
the squeezers or the 50:50 beam-splitters. 
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FIG. 13. (Color online) Steane Quantum Error Correction 
using the method of Glancy and Knill for shifts e~ lup which 
are detetected by the stabilizer check S q . The circuit only 
uses a single beam-splitter, a homodyne measurement of q 
and much stronger squeezing followed by a different corrective 
displacement D(fGK{q)) which is explicitly given in El- 


Appendix C: Choice of Feedback Phases 


Let us assume that the feedback phases tpx ,..., p m -i 
have been fixed and the bit string x[m — 1] = 
{xi... x m -{\ has been generated. How do we choose 
the next phase <p m in an optimal way? Using con¬ 
ditional probability distributions (i.e. the identities 
P(A\BC) = P(C\AB)P(A\B)/P(C\B) and P(C\AB) = 
P(CB\A)/P(B\A)), the probability distribution for 9 af¬ 
ter obtaining the measurement result a;[ to] = { x m , x[m — 
1]} equals: 


P tp (8\x[m]) 


P(9\x[m - 1 ])P Vm (x m | 6,x[m - 1]) 
P v (x m \x[m - 1]) 
P{0\x[m- l])P v (x[m\\9) 
P<p{x m \x[m - 1 })P(x[m - 1]|0) 
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Here we have explicitly kept the dependence of probabil¬ 
ity distributions on the phase ip which is chosen for the 
mth qubit measurement as this is the quantity that one 
wants to optimize over. This phase ip is chosen as the one 
which maximizes the sharpness of the a posterior proba¬ 
bility distribution P v (9\x[m\) averaged over possible out¬ 
comes x m (given the previous outcomes and choices for 
phases), or using Eq. (1T51) . 


<Pm = argmax ^ P v (x m \x[m - l])S'[P v ( 6 '|x[m])] 


z m = 0,1 


= argmax 


Xm — 0,1 


dde 


i8 P(Q\x[m - l])P v (x[m}\8) 


P(x[m — 1]|0) 


(C2) 


surement results of each round are independent, i.e. 

m m / f) 4 - \ 

P v=v > m {.x[m]\9) = Y[P v>i (x i \9) = ]Jcos 2 y-^r- +Xi \) ' 

(C3) 

If we assume no prior knowledge on 9 one can show 
that P(9\x[m — 1]) oc P{x[m — 1]|0) where the propor¬ 
tionality constant c is independent of 9 and ip m (but 
may depend on previous phases p m -i etc.). The ar¬ 
gument is as follows. For the first measurement we 
write P Vl ( 6 »|cci) = P lpi (xi\9)P(8)/P(x 1 ) oc P lfil (x 1 \9) 
if we assume that P(9) is a flat distribution (we have 
no prior knowledge about 9) and P(x i) = 1/2. Using 
Eq. m we see that this then holds for for all mea¬ 
surement outcomes by induction, i.e. P(9\x[m — 1]) oc 
P <p m ( x [ m ~ !]I 6 *)/ P Vm -i( x m-i\x[m - 2]). 

This implies that the optimization in Eq. (IC2I) is equiv¬ 
alent to 


ip m = arg max ^ 

v x m =0,l 


d9e w Pcp(x[m] |0) 


The expression for P v (x[m]\8) is simple as the mea- which can also be rewritten as 

I 


ip m = argmax 



dde ie 


cos 2 



Appendix D: Direct Controlled-Displacement 


We model the effect of applying a microwave drive to 
the cavity mode without including cavity decay. The rea¬ 
son for not including cavity decay is that the parameters 
of the protocol should be chosen to work without cavity 
decay such that the loss of photons from the cavity acts 
as a source of (shift) errors during this ideal protocol. 

The drive field is represented as a term Hdrive(i) = 
A (t)(a + as*) where A (t) = cos(a >dt) + H y (i) sin (ojdt) 

with ojd the drive frequency and fl y (t) the envelope 

of the pulse. The drive frequency ojd is chosen as u>d = 
oj r + x- The total Hamiltonian H to t{t ) = H e g + Hdrive(i) 
(in the lab frame) with H e ff in Eq. (12dl) equals 

= (uj r I -x(t) z )a^a + A(t)(a + a t ) - ( D1 ) 


In order to analyze the effect of the dynamics due 
to H to t{t), we use the following useful tool/fact. For a 
Hamiltonian of the form H(t) = uja'a + A (t)(a + at), the 
unitary time evolution can be written as a product of an 
overall rotation, an overall displacement and a phaseshift, 
i.e. 


p e ~iSo dt ' H(t‘) = /2( w r)D(7)exp(*$), (D2) 


-i / 0 T dt'\{t')e~ iut ' 






(C4) 


/ q T dt J) T dt'X(t)X(t') sin(w(t' — t)). We have 

R(ioT) = exp(— iujTa^a). This equality can be de¬ 
rived using the Suzuki-Trotter decomposition: 


n ( 

Te~^o TT I e ^A(t„)(a+at 

n—>oo ■*- 1 ^ s/ 

i=! V D t . 

= lim RZ(R- n D tn RZ)...(R~ 2 D t2 Rl)(R- 1 D tl R n ) 

n—> oo 

= R(coT)D ( lim V - iX ^ T e -i«Ti 

\ n —>oo ^' 77, 

V j=i 

where T is determined using D(a)D(f3 ) = 

exp(iIm(a/3*))P(a + /?) and taking the limit n —> oo. 

We can apply this equality to the cavity-transmon dy¬ 
namics U(0,T) = Texp (—i / 0 T dt'H tot (t')) as it is diago¬ 
nal in the (|0), |1)} basis: 

V(0 ,T) = R(( Ur - x)T)P(7-)e i t 3 '-+“« T / 2 t |0) (0| 

+ R((u r + x)T)D( 7+ )e i (*+-“* T/2 ) |1) (1|, 

where Ti = Jg dt J) T dt'X(t)X(t') sin((w r ± \)(t' — t)) and 
7 ± = —i / q T dt , A(t , )e _l t aJr±x ) t ) the Fourier transform of 




with 7 


and 
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A (t) at frequencies uj r ± We note that the relative 
rotation when qubit is in state |1) versus |0) is R{2\T). 
Thus without constraining \T, the interaction between 
qubit and cavity field realizes a controlled-rotation in 
addition to the desired controlled-displacement transfor¬ 
mation D( 7 ±). We can understand the undesired ef¬ 
fect of the relative rotation in the sketch in Fig. [T?] in 
which a controlled-displacement gate (and a some rela¬ 
tive controlled-rotation) is applied on, say, a cat state. 

If we demand that T is chosen such that \T = nk, the 
relative rotation acts trivially R(2\T) = I. One thus has 
to choose T = n/x which for x/27t = 2.5 MHz (which is 
the same as x/2-7t in the experiment in 0), equals 200 
nsec. 



FIG. 14. (Color online) Phase space sketches: the size of the 
blobs represents amplitude not quadrature uncertainty, (a) 
A cat state |a) + |—a) with a = y/n/2 (blue dots) gets dis¬ 
placed by ±\Jn/2 depending on a qubit state |0) or |1) along 
the 5 R(q) line producing the green dots. Assuming construc¬ 
tive interference of phases, the amplitude for |vac) is double 
that of the outlying states. Repetition of this scheme with 
constructive interference produces a binomial distribution of 
coherent states according to a Pascal triangle. In the phase 
estimation protocol, the input state is not a cat state but the 
squeezed vacuum (along the K(a)-direction) and the proto¬ 
col generates superpositions of displaced squeezed vacua on a 
line. In (b) it is shown how an additional phase space rotation 
which is different for |0) as for |1) produces a superposition 
of three states (red dots) which are no longer on a line, hence 
will not represent a code state. When the condition 2 \T = 2-k 
is met, this relative unwanted rotation vanishes. 


We can write down the expressions for the displace¬ 
ments D{ 7 ±). For a square displacement pulse (with 
fictitious zero rise time) which is turned on for the entire 
interval T, the displacements D{ 7 ±) are given by 


7+ = + S+: 

7- = , r T sincKr) 

+ H^ + ^) r e ixT sinc(xr) , 


(D3) 


with sinc(:r) = sin(a:)/a; and a displacement error c)_|_ = 

+ x)T)- We wish that 
7 _ ~ 0 (cavity is off-resonance) and note that this de¬ 
pends, in principle, on sinc(xT) and sinc(w r T). One 
can alway bound |sinc(a;)| < l/|x|, which shows that the 


dominant term in the expression for the off-resonant dis¬ 
placement 7 _ scales in strength like l s ^+^WI r . However, 
for our chosen pulse length T = 7r/x, sinc(xT) = 0 and 
this analysis is too pessimistic. 

For ui r /2TT = 10 GHz (for example) one has = 
^ = 2 x 10 3 so that |sinc(w r T)| < 0.5 x 10 -3 . For 
the well-justified approximation ui r x> l^+l ~ |7-| ~ 
| 7 +|sinc(u; r T) < | 7 + | x 0.5 x 10 -3 . This means that the 
off-resonant displacement is a factor 10 3 less than the 
on-resonant displacement. 

For the S p measurement, one can then choose Q x = 0 
and 7 + = y/2n = —%— (assuming 7 _ r; 0). Similarly, 
for the S q measurement one can choose Cl y = 0 etc. The 
value of the parameter Q y (or fl y ) depends on the exter¬ 
nal coupling of the cavity with the drive-line as well as the 
applied microwave power, e.g. Vt y oc y/n e xtFt where F t is 
the photon flux per unit time. For the square pulse this 
photon flux per unit time is assumed to be constant dur¬ 
ing the interval T. Thus P oc UdTFt = uidfipu i se is the to¬ 
tal input power and n pu ise the total number of photons in 
the displacement pulse. We thus find that one can enact 
the controlled-displacement gate when 2-7T a: K ex tn pu i se T. 
Given that T is fixed, depending on X; and ftext is re¬ 
quired to be small in order to have a high-Q cavity, this 
shows that by increasing n pu i se or the total input power 
one can always achieve the desired displacement. 

It is important to note that if the condition T = n/x 
is not accurately fulfilled, and one uses the pessimistic 
bound |sinc(:r)| < d, the relative difference in dis¬ 
placement equals | 7 _|/| 7 + | ~ = d. Also, it is 

not required that 7 _ « 0 as long as the value of the 
relative displacement (for, say, S p ) is large enough, 
i.e. 7 +^ 7 _ = y/2w and the values for 7 ± are known 
(calibrated). This means that we do not need to fullfill 
the condition T pu i se = T = n/x exactly. What happens 
when T pu i se < T and we also include the finite rise-time 
of the pulse? The expressions derived above will remain 
valid, but lead to a different y± (and 'F±) which could 
be fine-tuned using the pulse-shape Q, x {t),D, y {t). These 
values of 7 ± can also be calibrated by testing the 
controlled-displacement gate on the vacuum state. 

Let us comment on the phases which affect the 
transmon qubit and thus need to be taken into account 
if the qubit undergoes further rotations. Obviously, only 
the relative phase 'F + — will affect the qubit state. If 
this relative phase remains fixed for every round of the 
S p phase estimation measurement (and possibly different 
but fixed for the S q measurement), it will not affect the 
accuracy of the prepared code states, even if this relative 
phase is unknown. It essentially means that whenever 
one does controlled-<S , p in the protocol, one implements 
controlled- S p e^ + ~'^~' > instead. If the phase difference 
is not fixed, but known (so it represents a systematic 
error), it can be corrected with a rotation of the qubit. 
A fluctuating phase difference T + — will however lead 
to a noisier estimate in the phase estimation algorithm, 
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hence shift errors in the resulting code state. 
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FIG. 15. Probability P to prepare a code state with an effective error rate P^J r % using M ancilla qubits (rounds) and an 
infinitely squeezed vacuum state as initial state, binned to 0.2%. Errors due to finite squeezing are uncorrelated to the effective 
error rate shown here and can be taken into account using Fig. [8] The wide bin Pmoi% > 10% contains all measurement 
outcomes with effective error rate > 10%. Left: non-adaptive phase estimation of S p . Right: adaptive phase estimation of S p . 























































